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Abstract 

We consider a nonlinear evolution equation recently proposed to describe the 
small-x hadronic physics in the regime of very high gluon density. This is a func- 
tional Fokker-Planck equation in terms of a classical random color source, which 
represents the color charge density of the partons with large x. In the saturation 
regime of interest, the coefficients of this equation must be known to all orders in 
the source strength. In this first paper of a series of two, we carefully derive the 
evolution equation, via a matching between classical and quantum correlations, and 
set up the framework for the exact background source calculation of its coefficients. 
We address and clarify many of the subtleties and ambiguities which have plagued 
past attempts at an explicit construction of this equation. We also introduce the 
physical interpretation of the saturation regime at small Color Glass Con- 

densate. In the second paper we shall evaluate the expressions derived here, and 
compare them to known results in various limits. 
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1 Introduction 



There has been much activity in the last few years in an attempt to understand the physics 
of nuclear and hadronic processes in the regime where Bjorken's x becomes very small 



(see also ||2^ for a pedagogical presentation and more references). The remarkable 
feature about this regime is that the gluon density in the hadron wavefunction becomes so 
high that perturbation theory breaks down even for a small coupling constant, by strong 
non-linear effects. 

The important new phenomenon which is expected in these conditions is parton sat- 
uration |29[| , which is the fact that the density of partons per unit phase space 



{dN/(Pp± dyd'^x±) cannot grow forever as x becomes arbitrarilly small [here, y = ln(l/x) 
is the parton rapidity]. If this growth were to occur, then the cross section for deep inelas- 
tic scattering at fixed would grow unacceptably large and eventually violate unitarity 
bounds. One rather expects the cross section to approach some asymptotic value at high 
energy, that is, to saturate. It is also intuitively obvious that this should happen, since 
if the phase space density becomes too large, repulsive interactions are generated among 
the gluons, and eventually it will be energetically unfavorable to increase the density 
anyfurther. 

As X becomes smaller and smaller, the gluon density increases faster, and is the driving 
force towards saturation. This is the reason why, in what follows, we shall concentrate 
exclusively on the gluon dynamics. Saturation is expected at a gluon phase-space density 
of order l/og ^7\, 0, a value typical of condensates. This estimate, together with the 
fact that gluons are massless bosons, leads naturally to the expectation that the saturated 
gluons form a new form of matter, which is a Bose condensate. Since gluons carry color 
which is a gauge- dependent quantity, any gauge-invariant formulation will neccessarily 
involve an average over all colors, to restore the invariance. This averaging procedure 
bears a formal resemblance to the averaging over background fields done for spin glasses 



30|| . The new matter will be therefore called the Color Glass Condensate (CGC). This 
matter is universal in that it should describe the high gluon density part of any hadron 
wavefunction at sufficiently small x. 

The CGC picture holds in a frame in which the hadron propagates at the speed of 
light and, by Lorentz contraction, appears as an infinitesimally thin two-dimensional sheet 
located at the light-cone. The formalism supporting this picture is in terms of a classical 
effective theory valid for some range of x, or longitudinal moment umQ p"*". In this theory, 
originally proposed by McLerran and Venugopalan to describe the gluon distribution in 
large nuclei the valence quarks of the hadron (and, more generally, the "fast" partons, 
i.e., the partons which carry a large fraction of the hadron longitudinal momentum) 
are treated as sources for a classical color field which represents the small-x, or "soft", 
gluons in the hadron wavefunction. The classical approximation is appropriate since the 
saturated gluons have large occupation numbers ~ l/a^, and can thus be described by 

^See Sect. 2.1 for the definition of light-cone coordinates and momenta. 
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classical color fields with large amplitudes ~ 1/g. 

When going to smaller values of x, or lower longitudinal momenta p"*", one must 
integrate out the degrees of freedom associated with this change of scale ^ . This is done 
in perturbation theory for the quantum gluons (to leading logarithmic accuracy), but to 
all orders in the strong background fields which represent the condensate. This entails 
a change in the effective theory, which is governed by a functional, non-linear, evolution 
equation originally derived by Jalilian-Marian, Kovner, Leonidov and Weigert (JKLW) 



10, 11 1 . This is a functional Fokker-Planck equation for the statistical weight function 
associated with the random variable pa{x) (= the effective color charge density of the fast 
gluons) . 

One of the consequences of the JKLW renormalization group equation it that it pre- 
dicts the evolution of the gluon distribution function. In the low density, or weak-field 
limit, this equation linearizes and reduces to the BFKL equation as shown in Ref. 



TO]. But in the saturation regime, where A* ~ l/(yf, the JKLW equation is fully non- 
linear, and will be discussed in great detail in this and the accompanying paper (to be 
referred to as "Paper I" and "Paper II" in what follows). In particular, in Paper II, we 
will show that this leads to the Balitsky-Kovchegov equation in the large- A^^c limit (with 
Nc the number of colors). There is currently a discrepancy between the computations of 



Refs. [|12|, [T^, |T9|, ^ and those by Balitsky g] and Kovchegov ||2T| . The additional terms 
argued in do actually not show up in our computation [p2| . 

In the previous arguments, we have implicitly assumed the QCD coupling constant 
as to be small, which can be justified as follows: Associated with saturation, there is a 
characteristic transverse momentum scale Qs (a function of x), which is proportional to 
the density of gluons per unit area: 

where R is the radius of the hadronQ. Qs is the typical momentum scale of the saturated 
gluons, and increases with 1/x (see, e.g., Refs. [0, ^ |^). Thus, for sufficiently small 
^! Ql ^ ^QCDy the coupling "constant" asiQ"^) is weak indeed. (In any explicit 
calculation, the condition ^ ^qcd must be verified a posteriori, for consistency.) 

To summarize, the system under consideration is a weakly coupled high density system 
described by classical gluon fields. This is still non-perturbative, in the sense that the 
associated non-linear effects cannot be expanded in perturbation theory. But the non- 
perturbative aspects can now be studied in the simpler setting of a classical theory, via a 
combination of analytical and numerical methods (see, e.g., Refs. |jl|, |], ^ ^ ^ ^). In 
particular, in this context, it has been demonstrated [0, |^ that saturation occurs indeed, 
via non-linear dynamics. 

®This can certainly be defined for a nucleus, but even for a proton it is a meaningful quantity since, as 
we shall see, the typical transverse momentum at small x becomes large, so that one can properly resolve 
the transverse extent of the hadron. 
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The purpose of this paper and its sequel is twofold: first, to carefully rederive the 
JKLW evolution equation, by matching classical and quantum correlations; second, to 
render this equation precise and fully explicit, by computing its coefficients to all orders 
in the color field which describes the Color Glass Condensate. This requires an exact back- 
ground field calculation which will be prepared by the analysis in Paper I, and presented 
in detail in Paper II. 

Many of the results we shall present (especially in Paper I) are known in one form 
or another, and can be found in the literature. We shall nevertheless attempt to be 
comprehensive and pedagogical in our present treatment, in order to clarify several subtle 
points which have been often overlooked in the previous literature. To be specific, let us 
spell out here the new features of our analysis, to be presented below in this paper: 



The effective theory for the soft fields is stochastic Yang-Mills theory with a random 
color source localized near the light-cone (x^ ~ 0) [see Sect. 2.1]. In previous work, 
this source has been generally taken as a delta function at x~ = 0. Here, we find 
that in order to do a sensible computation, it is necessary at intermediate steps to 
spread out the source in x~ , that is, to take its longitudinal structure into account. 

For our 'spread out' source, we describe in detail the structure of the classical 
solution for the background field. We do this by relating the solution in the light 
cone gauge = 0) to the solution in the covariant gauge {dfj^A^^ = 0). We thus 
obtain an explicit expression for the background field in terms of the source in the 
covariant gauge [see Sect. 2.3]. This indicates that the source which is generated 
by quantum evolution must be evaluated in the covariant gauge as well. 

The quantum evolution of the effective theory is obtained by integrating out quan- 
tum fluctuations in a rapidity interval Ay = ln(l/6) (with 6 ^ 1) to leading order 
in ashi{l/b), but to all orders in the strong background flelds [cf. Sect. 3.2.1]. 
The quantum correlations induced in this way are then transferred to the stochastic 
classical theory, via a renormalization group equation of the Fokker-Planck type [cf. 
Sect. 3.2.2]. To prove the equivalence between the classical and quantum theories, 
we explicitly compute the induced fleld correlators {A^A^ . . .) in both theories and 
show that they are the same, to the order of interest [cf. Sect. 3.2]. 

The original JKLW renormalization group equation is written in terms of the color 



source in the light-cone gauge [TT|. When reexpressing this equation in terms of the 
source in the covariant gauge, we flnd that we must subtract out a term related 
to classical vacuum polarization, so that the classical and quantum theory properly 
match [see Sect. 3.3]. 

The non-linear evolution equation involves functional differentiations with respect 
to the source strength. With sources spread out and with explicit results from the 
covariant gauge, we have an explicit and transparent deflnition of these functional 
differentiations [see Sect. 3.2.2]. 
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The background fields are independent of tlie liglit-cone time , but inliomoge- 
neous, and even singular, in x~ . To cope with that, we find that it is more conve- 
nient to integrate the quantum fluctuations in layers of p~ (the light-cone energy), 
rather than of p"*" [cf. Sect. 3.4]. Accordingly, we have no ambiguity associated with 
possible singularities at p~ = 0. 

The gauge-invariant action describing the coupling of the quantum gluons to the 
classical color source is non-local in time UTO). Thus, the proper formulation of the 



quantum theory is along a Schwinger-Keldysh contour in the complex time plane 
[cf. Sect. 4]. However, in the approximations of interest, and given the speciflc 
nature of the background, the contour structure turns out not to be essential, so 



one can restrict oneself to the previous formulations in real time [10, 11 



We carefully flx the gauge in the quantum calculations. In the light-cone gauge, 
the gluon propagator has singularities at p"*" = associated with the residual gauge 
freedom under x~ -independent gauge transformations. We resolve these ambiguities 



by using a retarded ie prescription [cf. eq. (|3.84|) and Sect. 6.3], which is chosen 
for consistency with the boundary conditions imposed on the classical background 
field. We could have equally well used an advanced prescription. However, other 
gauge conditions, such as Leibbrandt-Mandelstam or principal value prescriptions, 
are for a variety of technical reasons shown to be unacceptable ^2 . 



• The choice of a gauge prescription has the interesting consequence to affect the 
spatial distribution of the source in x~ , and thus to influence the way we visualize 
the generation of the source via quantum evolution. With our retarded prescription, 
the source has support only at positive x~ [see Sect. 5.1]. 

• We verify explicitly that the obtention of the BFKL equation as the weak-fleld limit 
of the general renormalization group equation is independent of the gauge-flxing 
prescription [cf. Sect. 5]. 

Once these various points are properly understood, the background fleld calculation 
in Paper II becomes unambiguous and well deflned. As a result of this calculation, the 
coefficients in the JKLW equation will be obtained explicitly in terms of the color source 
in the covariant gauge |P2[. With such explicit expressions in hand, we hope that the 



functional evolution equation can be solved, via analytic or numerical methods, or both, 
at least in particular limits. In fact, as we shall see in Paper II, in the large Nc limit our 



results become equivalent to the equation proposed by Balitsky [ffl and Kovchegov El 



for which an analytic solution has been recently given [23 



The ultimate goal of this and related analyses is to derive the properties of hadronic 
physics at arbitrarily high energies from first principles, that is, from QCD. We conjecture 
that there will be universal behaviour for all hadrons in the high energy limit. The work 
we present here is a step in proving this, and in developing a formalism which allows for 
explicit calculations in the saturation regime. 
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The present paper is organized as follows: In Section 2, we discuss the solution to the 
classical Yang- Mills equations with a color source which define the Mclerran-Venugopalan 
model. In Section 3, we derive the non-linear evolution equation by matching classi- 
cal and quantum correlations. In Section 4, we formulate the quantized version of the 
Mclerran-Venugopalan model (on a complex-time contour), and derive Feynman rules for 
the coefficients of the renormalization group equation. In Section 5, we use the weak 
field, or BFKL, limit of the general evolution equations to test properties like the longi- 
tudinal distribution of the induced source, the equivalence between the renormalization 
group analysis in p"*" and in p~ , and the sensitivity to the ie prescription in the LC-gauge 
propagator. In Section 6, we exphcitly construct the gluon propagator in the presence of 
the background field describing the Color Glass Condensate. This propagator is the main 
ingredient of the quantum calculations to be presented in Paper II p2| . 



2 The classical approximation 

In this section, we shall introduce the McLerran-Venugopalan model [0] for the gluon 
distribution at small x. This is classical Yang-Mills theory with a random color source, 
which is the effective color charge of the partons having longitudinal momenta larger than 
the scale of interest. This simple model, which is motivated by the separation of scales 
in the infinite momentum frame (see Sect. 2.1), offers an explicit scenario for saturation 
0, |] (see Sect. 2.4) which supports the physical picture of the Color Glass Condensate. 
The MV model will be further substantiated by the analysis in Sect. 3, where we shall 
see that this model is consistent with the quantum evolution towards small x H, [10|, ITT 



2.1 The McLerran-Venugopalan model 

We consider a hadron in the infinite momentum frame (IMF), and in the light-cone (LC) 
gauge A'^ = 0, where the parton model is conventionally formulated (this gauge choice will 
be further motivated in Sect. 2.4). The hadron four-momentum reads P'^ = (P, 0,0,P) 
(we neglect the hadron mass), or, in LC vector notations, which we shall use systematically 
from now o 

= (p+,p-, P_l) = (v/2P, 0,0_l). The fast partons, that is, the hadron 
constituents which, like the valence quarks, carry a relatively large fraction of the hadron 
longitudinal momentum P"^, move as almost free particles with momenta collinear with 
P'^, and act as sources for soft gluons, i.e., gluons with small longitudinal momenta <^ 
P+, which are probed in deep inelastic scattering (DIS). 

For what follows, it is useful to have a sharp distinction between "fast" and "soft" 

^These are as defined as follows: for an arbitrary 4- vector v^, we write ~ {v^,v^,\-±), with 
w+ = (l/\/2)(t>" + v'^), = {l/\/2){v° — v^), and = {v^,v^). The dot product reads: p ■ x — 
p^x~ + p~ x'^ — p±_ ■ x± . p^ and p+ are, respectively, the LC energy and longitudinal momentum; 
correspondingly, x~^ and x^ are the LC time and longitudinal coordinate. 
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modes in the hadron wavefunction, according to their p+-momentum: to this aim we 
introduce an intermediate scale and define fast (soft) modes to have p"*" > A"*" (re- 
spectively, p"*" < A"*"). We choose A"*" = Xo-P"*", with xq not too small; in fact, the main 
constraint on xq is to be much larger than the Bjorken parameter x of the external probe 
that the hadron is interacting with. 

For instance, in DIS, this probe is the virtual photon 7*, with virtuality = —q'^q^- 
The Bjorken x variable is defined as x = Q^/2P ■ q and, by kinematics, it coincides with 
the longitudinal momentum fraction p'^ /P^ of the struck parton. In what follows, we shall 
always assume that x <^ 1. Rather than a virtual photon, which couples to the charged 
constituents of the hadron (quarks and antiquarks), it is more convenient to consider a 
DIS initiated by the "current" j = —\F^''F^^ which couples directly to gluons 0. Indeed, 
we are mainly interested here in the gluon distribution, which is the dominant component 
of the hadron wavefunction when x <^ 1. 

In the MV model, one assumes that the soft gluons can be described as the classical 
color fields radiated by the fast partons, themselves represented by a random color 
source p^. As we shall demonstrate in Sect. 3, this assumption is indeed justified (at least 
to lowest order in a^) provided the MV model is used as an effective theory valid for some 
range of soft momenta. Then, the fast degrees of freedom are treated as quantum, but 
they are integrated out (perturbatively) in the construction of the effective theory, while 
the soft modes are non-perturbative, because of their large occupation numbers, but can 
be treated as classical, for this same reason. Via this construction, to be detailed in Sects. 
3 and 4, the classical source pa will emerge as a natural way to describe the effects of the 
fast modes on the dynamics of the soft ones. In the original MV model this source 
has been simply postulated, and its properties have been inferred via an analysis of the 
separation of scales in the problem. Since this separation will play a major role in what 
follows, let us briefly discuss it here: 

The fast partons move along the z axis with large p+ momenta. They can emit, or 
absorb, soft gluons, but in a flrst approximation they do not deviate from their light-cone 
trajectories at z = t, or x~ = (no recoil, or eikonal approximation) . Thus, they generate 
a color current only in the -|- direction: = 5^^ Jf . 

As quantum flelds, the fast partons are delocalized in the x~ direction within a distance 
A~ ~ l/p"*"- However, when "seen" by the external probe or by the soft gluons (with 
momenta ^ P"^, and therefore a poor longitudinal resolution), they appear as sharply 
localized at the light cone, within a distance Ax^ ~ l/A"*". (This can be also understood 
as an effect of the Lorentz contraction). In some cases, to be carefully justifled later, the 
associated current J"*" can be represented as a 5 function at the hght cone: J"*" cx 5{x~). 
But in general, it turns out that the longitudinal extent of the source cannot be neglected: 
this will be important both for the classical calculations in this section, and for the 
quantum analysis to be performed in the subsequent sections, and in Paper II. 

Consider also the relevant LC-time scales: for on-shell excitations, 2p^p~' = p\, which 
implies that softer gluons have larger energies p^ ~ Q\lv^ 1 and therefore shorther life- 
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times Ax^ ~ ~ ~ x. Over such a short time interval, the "fast" partons are only 
slowly varying (since they have smaller energies), so their dynamics is essentially frozen. 
Thus, the soft gluons, or the external current in DIS, can probe only the equal-time cor- 
relators of the fast partons. All such correlators can be generated by a static^ current 
{quenched approximation). 

To summarize, the soft color current due to the fast partons is expected with the 
following structure: 

J^X) =(5'^+Pa(x-,X^), d-pa = ^ = 0, SUppp, = {X-|0<X- < 1/A+}, (2.1) 

which will be confirmed by the quantum analysis in Sect. 3-5 and Paper II. Note the 
restriction of the support of p to positive values of x~ . We shall see later that the precise 
longitudinal structure of the color source depends upon the gauge- fixing prescription, i.e., 
the condition used to completely fix the residual gauge freedom in the LC gauge. [Recall 
that, even after imposing the LC gauge, one still has the possibility to perform x~- 
independent (or "residual") gauge transformations, which preserve the condition A'^ = 
0. In quantum calculations, fixing this residual gauge freedom amounts to chosing a 
prescription for the pole of the gluon propagator at p~^ = 0.] With the specific prescription 
that we shall adopt (cf. Sects. 3.4 and 6.3), the source p is located at positive x~ , as 
anticipated in eq. (^.Ij). 

The current ( |2.1| ) acts as a source for the Yang-Mills equations describing the soft 
gluon dynamics: 

= 5^+p^ix-,^^). (2.2) 

To have a gauge-invariant formulation, the source pa must be treated as a stochastic vari- 
able with zero expectation value. This is also consistent with the physical interpretation 
of Pa as the instantaneous color charge of the fast partons "seen" by the shortlived soft 
gluons, at some arbitrary time. The (spatial) correlators of the classical variable Pa{x), 
with X = (a;~,xx), are inherited from the (generally time-dependent) quantum corre- 
lations of the fast gluons. The precise scheme for transforming quantum into classical 
correlations will be explained in Sect. 3. Here, we shall simply summarize the ensuing 
classical correlations in some (not yet specified) weight function for pa, denoted as PVa[p]- 
This is gauge-invariant by assumption, and depends upon the separation scale A"*" since it 
is obtained by integrating out quantum fiuctuations with longitudinal momenta | > A"*" 
(cf. Sect. 3). 

To conclude, in the MV model, the small-x gluon correlation functions are obtained 
in two steps: (i) First, one solves the classical Yang-Mills equations ( |2.2| ) in the light-cone 

*We will see below that, strictly speaking, this time independence assumption will have to be relaxed 
for a general solution of the classical equations of motion: The source turns out to be only covariantly 
time independent. We will however always be able to find a classical solution of the equations of motion 
which has a truly time-independent source. 
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gauge A'^ = 0. The solution will be some non-linear functional of p, which we denote 
as ^*(x)[p]. (Indeed, as we shall see in Sect. 2.3, one can always construct a solution 
which has A~ = and is time-independent.) (ii) The correlation functions of interest are 
evaluated on this classical solution, and then averaged over p, with the weight function 
Wa[p] : 

{Al,{x+,x)Ai{x-',y)...)A = Jvp WA[p]A:{x)Ai{y)-- - . (2.3) 
The normalization here is such as 

J Vp Wa[p] = 1. (2.4) 

Note that only equal-time correlators of the transverse fields can be computed in this 
way; but these are precisely the correlators which matter for the calculation of the gluon 
distribution function, cf. Sect. 2.4. 

Note furthermore that the correlations obtained in this way depend upon the scale A+. 
As we shall argue in Sect. 3 below, the effective theory in eqs. ( p.2| )-( p^ is valid only 
at soft momenta of order (that is, at momenta p+ ^ A+, but not too far below A+). 
Indeed, if one goes to the much softer scale bA~^ with 6 <^ 1, then there are large radiative 
corrections, of order asln(l/6), which invalidate the classical approximation at the scale 
A"*". In order to compute correlations at the new scale bA~^, one must first construct the 
effective theory valid at this scale, by integrating out the quantum degrees of freedom with 
longitudinal momenta in the strip bA~^ < \p'^\ < A^. Once this is done, the dependence 
upon the intermediate scale A"*" goes away, as it should. This will be explained in detail 
in Sects. 3 and 4. 

Remarkably, the equations ( |2.2| )-(|2l3|) above are those for a glass (here, a color glass): 
There is a source, and the source is averaged over. This is entirely analogous to what 
is done for spin glasses when one averages over background magnetic fields [^. To 
argue that one also has a condensate, one has to compute the correlation function (|2.3|) . 
By using a Gaussian approximation for the weight function, one has found a saturation 
regime where the classical field has a typical strength of order 1/g [0, ^, |2^ (see also 
Sect. 2.4 below). This is the maximal occupation number for a classical field, since larger 
occupation numbers are blocked by repulsive interactions of the gluon field. For weak 
coupling, this occupation number is large, and the gluons can be thought of as in some 
condensate. We are therefore lead to conclude that the matter which describes the small 
X part of a hadron wavefunction is a Color Glass Condensate. 

2.2 The Abelian case, as a warm up 

Before attempting to solve the Yang-Mills equation ( [^.2p , it is instructive to consider its 
Abelian version in some detail. This reads: 

d^F"'' = (5^+p(f) (2.5) 
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with F^"" = d^'A" - d^Af" and x = (x", x^). 

We are interested in LC-gauge (^4+ = 0) solutions which vanish as p ^ 0. Since 
p is static (i.e., independent of x~^), we can look for solutions which are static as well: 
d~A^ = 0. Then, the components of eq. (|2.5| ) with p = — and fi = i imply that A~ = 
(so that F~~^ = = 0), and also F^^ = 0. Thus, only the transverse fields A^{x) are 
non-vanishing, and define a "pure gauge" in the two-dimensional transverse planeQ. We 
can thus write [throughout this paper, we shall systematically use calligraphic letters (like 
A^ and to denote solutions of the classical equations of motion] : 

A\x) = -d'toix) (2.6) 

with uj{x) satisfying the following equation (cf. eq. ( p.5|) with fi = + and = —d^A^): 

- Vi(9+c.) = p (2.7) 

Eq. (|2.7|) can be easily solved in momentum space, to give: 

A'(p)=iMp) = -J^e^, (2.8) 

where, however, one needs a prescription to handle the pole at = 0. For instance, with 
the following, retarded prescription (with e 0+): 

1 _ 1 



(2.9) 



one obtains a solution which vanishes at x ^ — oo (thus, the "retardation" property 
refers here to x~ , and not to time) : 



A\x-,xj_)= dy-d'a{y-,x±), (2.10) 

J —oo 

with a{x) satisfying 

- Vla(f) =p(x). (2.11) 
This yields a^p) = p{p)/p\ (with p = (p^,p±)), or, in coordinate space, 

a{x) = - J ln(^|x_L - y±|p) p(2;",y^), (2.12) 

where p is an infrared cutoff which must disappear from physical quantities. 

In general, different prescriptions for the pole in l/p~^ correspond to different boundary 
conditions in x~ , which nevertheless describe the same physical situation, since they lead 
to the same electric field JF*"*" = —d~^A^ = d^a . Alternatively, solutions corresponding to 
different such prescriptions are related by (residual) gauge transformations. 

Note finally that if the external source is localized around x' = (cf. eq. (|2.1| )), so is 
a{x), and therefore the associated electric field jF*+(a;) = d^a{x). 

^Of course, this is not a pure gauge in four dimensions, since = 7^ —d^uj. 
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2.3 The non-Abelian case 



We now return to the non-Abelian problem, and note first that, written as it stands, 
eq. is not really consistent: The property [D^, [D^, F''^\\ = requires the color 

current J'^ in the Yang-Mills equations to be covariantly conserved, [D^, J^] = 0, which 
for a current = 5^+ J+ amounts to: 

[D-,J+] = d-J+ - ig[A-,J+] = 0. (2.13) 

In general, however, this is not satisfied by the static current (|2.1| ) (since d~ p = 0, but the 
commutator [^~,p] can be non-zero). Rather, eq. ( p.l3| ) shows that, for a non-zero field 
A~ , the current Jj" can be static only up to a color precession. Namely, if we identify pa 
in eq. (Ol) with the current at some particular time Xq, 



J+(x) = p(f) at x^ = x^, (2.14) 
then eq. ( p.l3| ) can be easily integrated out to give the current at some later time x~^ > Xq : 

J+(a;+,f) = W{x^,x) p{x)W\x^,x). (2.15) 
We have introduced here the temporal Wilson line: 

W[A-]{x+,x) = Texpjic/^^ dz+yl~(z+, f)| , (2.16) 

with T denoting the time ordering of the color matrices in the exponential (that is, the 
matrix fields A~{x) = A~{x)T°- are ordered from right to left in increasing sequence of 
their x~^ arguments). 

The ensuing equations of motion for the background field: 

[D^,F''^']{x) = 6^'+W{x)p{x)W\x), (2.17) 

are generally complicated by the non-locality of the current in time. As in the Abelian 
case, however, it is consistent to look for solutions to eq. ( |2.17| ) which are staticQ and 
satisfy A~ = (in addition to the LC gauge condition A'^ = 0): 

A^ = A- = 0, A' = A\x-,x±). (2.18) 

This Ansatz is preserved by gauge transformations which are both x~- and ^"'"-independent, 
that is, by two-dimensional gauge transformations in the transverse plane. With this 
Ansatz, eq. (|2.17| ) reduces back to eq. ( |2.2| ), which explains the emphasis put on this 
equation in Sect. 2.1. In particular, for p = +, 



[A,F*+] = p{x), (2.19) 



^"^Of course, time-dependent solutions with non-z ero A „ can be also constructed, e.g., via ^"'"-dependent 
gauge transformations of the static solution in eq. ( 2. IS ). For what follows, it is sim ply more convenient 



to use the residual gauge freedom in order to choose the solution in the form ( p.l 
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while the ^ = i component, [Dj, F^^] = 0, imphes that is a pure gauge in the transverse 
plane (JF^'* = 0). We thus write: 

A\x-,x±) = ^U{x',x^)d'U\x-,x±) , (2.20) 

with the SU{N) group element U{x~,x±) imphcitly related to p by eq. ( |2.19| ). 

To make progress, it is useful to observe that the fields in eq. ( p. 201 ) can be gauge 
transformed to zero by the gauge transformation with matrix W{x): 

A^' — > A^" = U^Af'U + - U^d^U. (2.21) 

9 

This yields indeed = A^ = 0, together with 

i+ = 1 f/t (^d+u) . (2.22) 

Thus, there exists a gauge where the non-Abelian field has just one non-trivial component, 
A'^ = S'^'^A'^. This gauge will play an important role in what follows. It will be referred to 
as the covariant gauge, since d^A^ = for the fields in eqs. ( |2.21| )- (|2.22| ). In this gauge, 
the general Yang- Mills equations ( p.l?| ) reduce to a single linear equation, 

-Vli+(f) = p(x), (2.23) 

where 

p{x) = U^{x) p{x)U{x) (2.24) 

is the classical color source in the covariant gauge. (We preserve the notations without 
the tilde — e.g., p — for quantities in the LC gauge.) 

Eq. ( p.23| ) is formally the same as eq. ( |2.11| ) for a in QED, and we prefer to rename 
A~^{x) = a{x) in what follows. That is, a satisfies 

- Via(f) = p(f) , (2.25) 

and is a linear functional of the covariant gauge source p. Thus, the classical solution 
in the covariant gauge is trivial to obtain. However, physical quantities like the gluon 
distribution function rather involve the correlators of the fields A^ in the light-cone gaugeQ 
(cf. eq. ( P73| ) and Sect. 2.4). Still, the fact that the latter can be expressed as gauge 
rotations of the quasi-Abelian field A^ = will greatly simplify the calculations to 

follow. Specifically, eq. (|2.22|) can be inverted to give (with A'^ = a) 

{x~ , x±) = P exp lig J dz~ a{z~ , x±)\ , (2.26) 



^^Of course, physical quantities are gauge invariant, and could be computed in any gauge; in gen- 
eral gauges, however, they are not simply related to the Green's functions of the corresponding vector 
potentials; see also Sect. 2.4. 
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where the symbol P orders the color matrices a{x) from right to left, in increasing or 
decreasing order of their x~ arguments depending respectively on whether x' > Xq or 
x~ < Xq . The initial point yet arbitrary, and will be specified in a moment. 

(Different choices for Xq lead to LC-gauge solutions which are connected by residual, 
two-dimensional, gauge transformations.) 



Together, eqs. ( |2.2CI| ), ( |2.25| ) and (|2.26|) provide an explicit expression for the LC- 



gauge solution as a non-linear functional of p, to be denoted as A^[p\. If one rather 
tries to express A^ in terms of p, then the corresponding functional A^[p] would be more 
complicated and not explicitly known. Indeed, a is related to p via the following non-linear 
equation (cf. eqs. (g2^)-(g^)y. 

U[-Vla)u^ = p, (2.27) 

for which we have no explicit solution. 

For what follows, it will be essential to know the solution A^ explicitly. We thus 
propose to use the covariant gauge source p, rather than the LC-gauge source p, as the 
functional variable that the fields depend upon. This is possible because the measure and 
the weight function in eq. (|2.3| ) are gauge invariant, so that the final average over the 
classical source can be equally well expressed as a functional integral over p. In other 
terms, from now on we shall replace eq. (p.3|) with 



{A\x+, x)A^(x+, y) ■ ■ ■)a = / Vp Wa[p] AUp] Al[p] ■ ■ ■ . (2.28) 

Up to this point, the spatial distribution of the source in x~ has been completely 
arbitrary: the equations above hold for any function p{x~). For what follows, however, 
it is useful to recall that p has the localized structure shown in eq. (|2.1|) . As mentioned 
after this equation, and it will be demonstrated by the quantum analysis to follow, the 
precise longitudinal structure of p is sensitive to the conditions used to completely fix the 
gauge, in both the classical and quantum calculations. 



To fix the classical gauge, i.e., the gauge for the classical solution (|2.2CI| ), we shall 



adopt retarded boundary conditions in x~ : A^{x) — > for x~ — oo . That is, we choose 
Xq — i> — oo in eq. ( p.26|) . This amounts to a complete gauge fixing since this boundary 
condition would be violated by any ^"-independent gauge transformation. Remarkably, 
we shall see later that this classical gauge condition also fixes the gauge prescription to 
be used in the quantum calculations: indeed, for the above boundary condition to be 
consistent with the quantum evolution, one must adopt a retarded ie prescription for the 
axial pole in the gluon propagator (see Sect. 6.3). With these gauge-fixing prescriptions, 
the ensuing classical source has support only at positive x~ , with ^ x~ ^ l/A"*" (cf. 
Sect. 5.1 and Paper II). 

In what follows. A"*" will be most often the large longitudinal momentum scale in the 
problem: both the external current in DIS, and the quantum fluctuations to be considered 
in Sect. 3, will have momenta ^ A^. Because of their poor longitudinal resolution. 



12 



such fields are sensitive only to the gross features of the background fields over large 
distances ^ where we can write: 



U\x',x±) = Pexpji^ J dz- a{z",x±)j ^ e{x')V\x±) + e{-x'), (2.29) 

with 

V^{x^) = Pexpjz^ dz' a{z-,x^)Y (2.30) 



Together with eq. (|2.20 ), this gives a background field of the form: 



A'{x-,x^) ^ 0{x')-V{d'V^) = e{x-)A'^{x^). (2.31) 

The associated electric field strength is then effectively a (5-function: 

r^{x) = -d+A' ^ -6{x-)A'^ix±). (2.32) 

It is implicitly understood here that the various 6- and ^-functions are regularized over 
a distance Ax~ ~ 1/A+. 

2.4 Gluon distribution function and saturation 



According to eq. (|2.28| ), there are two ingredients needed in order to compute soft gluon 



correlations in the MV model: the solution A^^lp] to the classical equations of motion 
and the weight function PVa[p]. The classical solution has just been constructed, and is 
given by eqs. ( |2.20| ), ( p.25| ) and ( |2.26| ) (with x^ — oo). The construction of the weight 
function from quantum fluctuations will be the main objective of the remaining part of 
this paper together with Paper II. But before turning to the quantum dynamics, let us 
briefly recall the results of some previous calculations in the MV model 0, |^, which use 
a Gaussian weight function and show saturation explicitly. This will also give us the 
opportunity to discuss the gauge-invariant definition of the gluon distribution function. 

We first recall the usual definition of the gluon distribution function, which is written 



in the LC gauge |26 



Ai(x+, k*, kx)^' (1+, -k+, -kj_)), (2.33) 
with the brackets denoting an average over the hadron wavefunction. Eq. ( p.33| ) can be 



understood as follows ||2^: In light-cone quantization, and with k = (fc+,k_|_). 



^ A^(x+, k)J^^{x\ -^) = E E 4c(A:) a,.{k) = — (2.34) 
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is the Fock space gluon distribution function, i.e., the number of gluons per unit of volume 
in momentum space. (Here, a\^{k) and axc{k) are, respectively, creation and annihilation 
operators for gluons with momentum k, color c and transverse polarization A, and we use 
the normalization in Ref. |^.) Thus, eq. ( p. 331) simply counts all the gluons in the hadron 



wavefunction having longitudinal momentum k~^ = xP^, and transverse momentum up 
to Q. 

As emphasized above, eq. (|2.33| ) is meaningful only in the LC gauge ^4+ = 0; indeed, it 
is only in this gauge that the matrix element in the r.h.s. has a gauge-invariant meaning, 
as we discuss now: Note first that, in this gauge, the electric field F*"*" and the vector 
potential A^^ are hnearly related, F^'^{k) = ik^Al^{k), so that eq. (|2.33|) can be also written 



as a two-point Green's function of the (gauge-covariant) electric fields. After performing 
the integral over k^, one obtains (with k^ = xP^ from now on) : 

xG{x, Q') = lJ^, QiQ' - kl) {Ft{x\ k)F:^{x\ -k)). (2.35) 

This does not look gauge invariant as yet: in coordinate space (a; ■ = a;~A;+ — x_|_ ■ k_|_), 

F^+(x+,fc)F^+(x+,-fc) = jd'xj Sy e^(^^-^^)-^>^+(x+,f)F^+(x+,f) (2.36) 

involves the electric fields at different spatial points x and y. A manifestly gauge invariant 
operator can be constructed by inserting Wilson lines: 

Tr {F*+(f) f/^(f , y) F'+{y) U,{y, x)} , (2.37) 

where (with A = (A+, A_l)) 

U^{x,y) = Pexpji^ J dz-A{z)Y (2-38) 

and the temporal coordinates x~^ are omitted (they are the same for all the fields). In 
eq. ( |2.38| ), 7 is an arbitrary oriented path going from y to x. For any such a path, eq. ( |2.37] ) 
defines a gauge-invariant operator. 

In particular, let us choose a path made of the following three elements: two "vertical" 
pieces going along the x~ axis from {y~, y±) to {—00, y±), and, respectively, from {—00, x±) 
to {x~,x±), and an "horizontal" piece from {—oo,y±) to {—oo,x±). Along the vertical 
pieces, dz-A = dz~A~^, so these pieces do not matter in the LC gauge. Along the horizontal 
piece dz ■ A = dz± ■ A^(— 00, z±), and the path 7 is still arbitrary, since this can be any 
path joining y± to x± in the transverse plane at x~ —00. This arbitrariness disappears, 
however, for the classical solution constructed in Sect. 2.3: this is a two-dimensional pure 
gauge (cf. eq. ( p^.20| )), so the associated Wilson lines in the transverse plane are path- 
independent. Moreover, these Wilson lines become even trivial, U^{x±,y±)x-,y-^-oo = 1) 
once we choose the retarded prescription Xq —>■ —00 in eq. ( p.26| ). 
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Thus, within the classical approximation, for the particular class of paths mentioned 
above, and with retarded boundary conditions for the classical solution, the manifestly 
gauge-invariant operator product in eq. (2.37) coincides with the simpler operator which 



enters the standard definition of the gluon distribution, in eqs. (|2.35|) - (|2.36|) . Converserly, 
the latter can be given a gauge- invariant meaning under the conditions mentioned aboveQ. 

In the remaining part of this section, we shall be concerned only with this classical 
approximation, which is expected to work better at very small x, in particular, in the 
saturation regime. Then we can replace, in eq. ( p.35|) , F^'^ —* J^^'^, with JF*+ the time- 
independent classical electric field constructed in Sect. 2.3 (a functional of p): 

xG^,{x,Q') = IJ 0^e(g2-fcl)(|J-:+(^)p)^. (2.39) 

Here and from now on, the average is to be understood as an average over p, in the sense 
of eq. ( ^.2(j| ) where the scale A"*" is now chosen as A+ ~ xP+ (cf. the discussion after 
eqs. (pD-(|3)). With eq. (|3D for J^^+, eq. (pD reduces to: 



= R"f^2 J rfVe-^'^---(^-(0)^-(x^))^, (2.40) 

where R is the hadron radius, and we have assumed homogeneity in the transverse plane, 
for simplicity. 

Note that, if there were not for the x-dependence hidden in the weight function for p 
(that is, PVa[p] with ~ xP^), the r.h.s. of eq. (|2.40|) would be independent of x, and 



so would be also the quantity xGci{x,Q'^), in agreement with leading-order calculations 
in light-cone perturbation theory [^. Thus, in the MV model, the actual x-dependence 



of the gluon distribution function is encoded in the weight function, and is a consequence 
of the quantum evolution (which makes W[p] a function of A"*"; cf. Sect. 3). 

To make progress, a model for the weight function Wa[p] is required. As argued in 
Refs. 0, a simple approximation which should be valid for transverse wavelengths much 
smaller than the size of the hadron, is to write 



H--.|pl==exp|--yd^x|^|. (2.41) 

where ^\ is proportional to the total color charge density squared of the partons with 
p'^ > A"*". For instance, in a simple valence quark model, 

d(^)-/f^. (2.42) 



' c 



""^^Note also that, if we let ^ cx) in eq. ( 2.35 ), then the corresponding limit of the gluon distribution 



is gauge invariant in full generality, since the unrestricted integral over k± sets x± = y± 
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where n{x ,x±) is the quark number density in the hadron, normahzed such as (for a 
nucleon) 

dx~ J d'^x±n{x~,x±) = Nc. (2.43) 

With this approximation for PVa[p], let us first compute the gluon distribution function 
( |2.39 ) in the linear response approximation, that is, for a source p which is weak enough 



so that one can linearize the results in Sect. 2.3 (this is the long as X is not too 

small). One thus obtains JF^-' ^ i{y /k'jjpa , cf. eq. (p.8|), and therefore: 

xG.(x.O^) i /|^ ^"^lr''' (lP-Wr)v (2-44) 

(Note that there is no difference between p and p in this linear approximation.) The 
average over p is now easily performed by using eqs. ( p.41| )-( |2.43| ), which give (with 
as = gy^TT and Cp = (N^ - l)/2N,) : 

xGM,Q') ^^^^ln%, (2.45) 



where p is an IR cutoff as in eq. ( p.l2|) . This result is as expected from LC perturbation 
theory |l26 |. 



But in the present formalism, one can actually do better: with the Gaussian weight 
function (|2.41| ) and the non-linear classical solution in Sect. 2.3, one has been able to 



evaluate the fully non-linear expectation value in eq. ( |2.40| ), with the following result 

AZ(^.)} = ^ ;j . (2.46) 

where Qs oc a^^A is the saturation momentum, and is expected to increase when A+ (or 
x) decreases. (The above equation is valid only for x± <^ 1/Aqcd- The logarithmic 
dependence upon Aqcd comes from properly regulating the infrared transverse momenta, 
and is associated with the over all color neutrality on scale sizes of order 1/Aqcd 0-) 

The remarkable feature about eq. ( |2.46D is that it displays saturation via non-linear 
effects. Namely, the LC gauge potential never becomes larger than 

As anticipated in the Introduction, gluon saturation requires fields as strong as A* ~ l/(yf, 
which supports the physical picture of a condensate. Together with eq. ( |2.2| ), this implies 
that p ~ 1/(7 at saturation. 
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This interpretation can be made sharper by going to momentum space. As obvious 
from eq. ( p.40|) , 



N{k^) ^ ' J> ^ / rfVe-^^---a-(0)^-(xx)), (2.48) 



d?k±^ d?'h±_ 

is the gluon density per unit of x and per unit of transverse phase-space. By using this 
and eq. ( |2.46[ ), one obtains |^ 

N{k^) oc as{Ql/kl) for kl > Q], (2.49) 

which is the normal perturbative behavior, but 

lln^' 

as Ql 



N{k^) oc — In for kl < (2.50) 



which shows a much slower increase, i.e., saturation, at low momenta (with /c^ ^ ^qcd 
though). 

Note, however, that the above argument is not rigorous, since the local Gaussian 
form for W in eq. ( p.41| ) is only valid at sufficiently large transverse momentum scales 
so that the effects of high gluon density are small. It is therefore important to verify 
if saturation comes up similarly with a more realistic form for the weight function, as 
obtained after including quantum evolution towards small x. This would also determine 
the a;-dependence of the saturation scale. We now turn towards this quantum analysis. 



3 The non-linear evolution equation 

With this section, we begin the study of the quantum dynamics which allows one to 
actually derive the MV model as an effective theory valid for some range of longitudinal 
momenta. The central element in this construction is a functional renormalization group 
equation (RGE) which describes the flow of the weight function PVa[p] with the separation 
scale A+. By means of this equation, the soft correlations induced by the quantum fields 
with momenta p+ > A+ are transferred into classical correlations of the random source p. 
Thus, this equation governs the quantum evolution of the effective theory. 



In Refs. |jTO|, 0, this RGE has been obtained by studying the correlators {6p6p. . .) 
of the color charge 6p of the quantum fluctuations. However, physical quantities like the 
gluon distribution function are more directly related to field correlators, like {A^A^ . . .) 



(see, e.g., eq. ( |2.33|) ). Given the different ways that the fleld correlators develop out of 
the charge correlators in the classical and quantum settings, it is not a priori obvious 
that matching charge correlations should be enough to insure the equivalence between 
the two theories. Our approach here will be therefore different: the RGE will be derived 
by matching directly fleld correlators, to the accuracy of interest. This will lead us indeed 
to the JKLW evolution equation for VFa[p] when the source p is in the light-cone gauge. 
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In addition, we shall also establish the evolution equation for VFa[p] where the source 
p is in the covariant gauge. According to the discussion in Sect. 2.3, it is this latter 
representation which is more suitable for explicit calculations. 

3.1 A quantum extension of the MV model 

To study the quantum evolution of the MV model, one needs first a quantum extension of 
this theory. By definition, such an extension should allow for soft < A"*") quantum 
gluons in addition to the classical color fields generated by p. Only soft fluctuations 
are permitted since, by assumption, the fast > A+) degrees of freedom have been 

already integrated out in the construction of the effective theory at the scale A"*". 

A generic quantum extension of the MV model is defined by the following generating 
functional for soft gluon correlation functions: 



This is written fully in the LC gauge (in particular, p is the LC-gauge color source), and 
involves two functional integrals: a quantum path integral over the soft gluon fields : 



(with Zx[p] = Z/y[p,j = 0]), which defines quantum correlations at fixed p, and a classical 
average over p with weight function PVa[p]. The upper script 'A" on the quantum path 
integral is to recall the restriction to soft {\p'^\ < A"*") longitudinal momenta. Note that 
the separation between fast and soft degrees of freedom according to their longitudinal 
momenta has a gauge-invariant meaning (within the LC-gauge) since the residual gauge 
transformations, being independent of x~ , cannot change the momenta. 

All non-trivial information about the quantum dynamics is contained in the action 
S[A, p]. This is chosen such as to reproduce the classical equations of motion ( p. 171) in the 
saddle point approximation 6S/6A^ = (so that, at tree-level, the quantum theory ( |3.1| ) 
reduce to the classical MV model, as it should). A gauge-invariant action satisfying this 
property has been proposed in Refs. [T^, ITT], (see also Ref. ^i^^ t)e presented 

in Sect. 4 below. As we shall argue there, this original proposal requires some further 
refinement (in the form of a complex-time formulation), but this is not essential for the 
general discussion in the present section. In fact, for what follows, it is enough to know 
that S has two pieces, 5'[y4, p] = S'yA/[^] + Sw[A^ , p], with Sym the usual Yang-Mills 
action, and Sw a gauge-invariant generalization of the eikonal vertex / d^xpaA~. (See 
Sect. 4 for more details.) 

At tree-level, A'^ ^ A'^ = S^''Al_, with A^ the classical field generated by pa, cf. Sect. 
2.3. In general, the total gluon field in eqs. (p.l|) - (p.2|) can be decomposed into tree-level 
field plus quantum fluctuations: 



z[j] = jvp WM z-Ap\ I i^A>;^s{At) e^^t^'^i-/^-^ 



(3.1) 



(3.2) 



A^^ix) = A^^{x) + SA^^{x). 



(3.3) 
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The average field in the system is not just A'^, but involves also an induced piece {SA'^) 
coming from the polarization of the quantum gluons by the external source : 



(A^^ix)) = A^^ix) + ^ A^^x) + 6A^Ax). (3.4) 

This average field satisfies the following equation 

/ SS \ 



6A{i(x) 



0, (3.5) 



where the brackets refer to the average over quantum fiuctuations at fixed p. (Unless 
otherwise specified, we shall always use this convention in what follows.) 



From eq. ( |3.1| ), the (time-ordered) two-point function of the soft gluons is obtained as: 

« T..(.).,.) » - W.l, {tlA^M^f^)^ , ,3.) 

where we use double brackets to denote the average over both quantum fiuctuations and 
the random source. This double average should be contrasted with 

fVp WaIp] I'^VA A^(x)A-(y) e-^[-^-^] 

JVp WaIp] J^VA e^^I^'H ' ^ ■ ^ 

which would be the correct average if p and A^ were two sets of quantum fields which are 
coupled by the dynamics (e.g., fast and soft gluon fields). In reality, Pa{x) is a classical 
variable which has no dynamics by itself, and represents only in an average way the effects 
of the fast partons on the dynamics of the soft gluons. Then the appropriate averaging 



is that in eq. (3^), where one first performs a quantum average at fixed p, and then 
a classical statistical average over p. This refiects the physical fact that the changes of 
p happen on a time scale which is much larger than the time scale characterizing the 
dynamics of the soft gluons. This is very similar to the average over disorder performed 
in the context of amorphous materials (hke spin glasses) [Q, and supports the physical 
picture of the saturation regime as a Color Glass Condensate. 

Since the intermediate scale A is arbitrary, it must cancel out in any complete cal- 
culation of soft correlations. That is, the cutoff dependence of the quantum loops must 
cancel against the corresponding dependence of the classical weight function PVa[p]- This 
constraint can be formulated as a renormalization group equation for VFa[p] ||T1| , 0, to 
be constructed in the next subsection. 



3.2 From quantum to classical correlations 



In writing down the quantum effective theory in eq. (|3.1|) , we have assumed that the 
effects of the fast gluons can be reproduced — at least to some accuracy which, as yet. 
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has not been specified — by a classical color source p with a peculiar structure (time- 
independent, and localized near the light-cone), and some (undetermined) weight function 
PVa[p]- In this subsection, we show how to construct this effective theory, step by step, 
by integrating quantum fiuctuations in successive layers of p"*" (or p~). 



As usual with the BFKL evolution |31|, the most important quantum corrections at 
small X are those which are enhanced by large intervals of rapidity Ar = ln(l/a;) 3> 1. 
As we shall see in the explicit calculations in Sect. 5 and Paper II, the amplitudes with 
soft {k^ < A"^) external lines receive quantum corrections of order ln(A"*"/A;+) from 
integration over the fiuctuations with p"*" momenta in the strip fc"^ <^ ^ A+. (A"*" is 
the upper cutoff on the longitudinal momenta of the quantum fiuctuations, cf. eq. ( |3.1| ).) 
Clearly, this is a large correction as long as A"*" ^ fc"*". Converserly, if A"*" is small enough 
(close to fc"*"), the quantum effects are relatively small (ordinary perturbative corrections 
of order a^), and the classical approximation is justified, at least, to lowest order in a^. 

This argument suggests that one can use the classical MV model to study soft correla- 
tions provided all the fast partons with momenta larger than the scale of interest have 
been integrated out, and their effects included in the parameters of the classical theory. In 
this way, all the potentially large logarithms are resummed in the structure of the weight 
function H/a[p], with A"*" ~ . The resulting classical theory is then the correct effective 
theory at the scale fc"'", up to corrections of order (with no large logarithms). 

To construct the classical theory, we shall study the evolution of p and PVa[p] with 
decreasing A+. Specifically, we shall consider a sequence of two effective theories ("Theory 
I" and "Theory 11") defined as in eq. (|3.1|) , but with different separation scales: A"*" in 
the case of Theory I, and feA"*" for Theory II, with 6 <C 1, but such as Q;sln(l/6) < 1. 
(This allows us to treat corrections of order a^j ln(l/6) in perturbation theory.) Theory II 
differs from Theory I in that the quantum fluctuations with longitudinal momenta inside 
the strip 

5A+ < < A+, (3.8) 

have been integrated out, and the associated correlations have been incorporated at tree- 
level, within the new weight function W^p^. 

The difference AVT = PF^a — W^a will be obtained by matching calculations of gluon 
correlations at the scale k~^ ^ hK^ in the two theories. In Theory II, and to lowest order 
in as-, these correlations are found already at tree level, i.e., in the classical approximation. 
In Theory I, and to the same accuracy, they involve also the logarithmically enhanced 
quantum corrections due to fluctuations with momenta in the strip. In computing the 
latter, we shall work to leading order in ln(l/6) ["leading logarithmic accuracy" (LLA)], 
but to all orders in the background flelds and sources. Indeed, we are interested here in 
the regime of saturation, where ~ 1/ g and p ~ ^/ 9, so that the associated non- linear 
effects cannot be expanded in perturbation theory. 

The final outcome will be a functional, and non-linear, differential equation for VFt-[p] 
(with r = ln(l/6)) which describes the flow of the weight function with r [10, 11]. In 
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the weak field limit, this equation can be linearized, and shown to reduce to the BFKL 
equation |jTO| (see also Sect. 5). 



3.2.1 Quantum corrections in Theory I 

The classical effective theory is expected to reproduce equal-time correlators of the trans- 
verse fields with soft {k~^ ^ longitudinal momenta (cf. eq. (|2.3|) ). So, these are the 
correlations that we shall try to match between Theory I and Theory II. In turns out 
that, in order to establish the evolution equation for Wr[p], it is enough to consider the 
two-point function 

{Al{x+,k)Al,{x+,-k)), (3.9) 

where k = {k~^, k_|_) with k~^ ^ bA'^, and the brackets denote quantum expectation values 
at fixed p. This is independent of time because the background is static. When computing 
this function within Theory I, there are quantum corrections of order asln(l/6) to be 
identified in what follows. To simplify writing, it is convenient to work in the coordinate 
representation, and define (with A^. = v4^(a;)) 

Gix^y) ^ (Al(x+,x)Al(a;+,y)) = {{A' + 6A%{A' + 6A%) , (3.10) 

where it is understood that the external lines carry soft (^ 6A+) longitudinal momenta. 

The equal-time correlator (|3.9| ) or (|3.10|) represents the density of the soft gluons 
generated by p (cf. Sect. 2.4). [The vacuum, p-independent, contributions to this quantity 
are to be subtracted away] It involves a tree- level piece A\x)A\y) together with quantum 
corrections — the induced mean field dA^ = {SAD, and the induced density {6Al^6AD 
— which describe the polarization of the quantum fiuctuations in the presence of p. In 
general, these corrections receive contributions from all the fiuctuations with momenta 
\p'^\ < A"*". However, for matching purposes, we need only the respective contributions of 
the semi-fast gluons, by which we mean the fiuctuations with p"*" momenta inside the strip 
( |3.8[ ). Indeed, these are the only quantum effects from Theory I that should be found at 
the tree-level of Theory II. (The remaining quantum effects, due to the fluctuations with 
Ip^I < 6A^, will appear as radiative corrections also in Theory II.) 

From now on, we shall include in 6A^ or {6A^6A^) only the quantum effects generated 
by the semi-fast gluons. In fact, to LLA, the definition of the "semi-fast gluons" — 
i.e., the quantum fiuctuations that have to be integrated out in going from Theory I to 
Theory II — can be restricted even further: These are the nearly on-shell fiuctuations 
with longitudinal momenta deeply inside the strip, 6A+ -C \p'^\ <^ A"*", and energies 
A~ ^ Ip^I ^ ^ where 

A--^. (3.11) 
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and Q±_ is a typical transverse momentum. (To LLA, it makes no difference what is 
the precise value of Q^.) Indeed, as we shall see later, these are the modes which are 
responsible for the logarithmic enhancement of the radiative corrections. 

For more clarity, we introduce the notation a^(a;) for the semi- fast gluons, and preserve 
the notation 5A^{x) for the softer fields, with momenta \k'^\ ^ (The fields in 

eq. ( p.lO|) belong to the latter category.) That is, we rewrite the total gluon field as 



follows (compare to eq. (|37^ )) : 



A(^(x) = A^:,{x) + 5A^Sx) + <(x). (3.12) 

The relevant quantum effects arise from the interactions between the soft fields 5A^ and 
the semi-fast gluons in the background of the tree-level fields and sources and p. 
This yields, as we shall see, 5 A' ~ a,ln(l/6)^* and {5 A' 5 A') ~ ln(l/6)^*^*, so that 

g{x, y) = A\x)A\y) + A\x)6A\y) + 6A\x)A'{y) + (M'(x+, x)6A\x+, y)), (3.13) 

where the disconnected piece 6A^SA^ of {6A^6A^) can be discarded, since of higher order in 
as- In general, effects non-linear in 6A^ can be ignored since the induced fields, in contrast 
to the tree-level fields, are weak. Because the background fields A'' and p are static, so is 
the induced mean field 6A^ (and any other one-point function), while two-point functions 
like {6Al.6Ajj) depend only upon the time difference x~^ — ■ In particular, (a^) = 0, 
since the field has no static mode. 

In what follows, we shall express the quantum effects SA and {5A^5A^) in terms of 
correlation functions of the semi-fast gluons. In doing this, we shall perform simplifica- 
tions appropriate to LLA. First, the semi-fast gluons will be treated in the Gaussian, or 
mean- field, approximation, which means that we shall consider, at most, one- loop dia- 
grams (higher loops contribute to higher orders in as)- Then, we shall perform kinematical 
approximations which exploit the hierarchy of scales to retain only terms of leading or- 
der in asln(l/6). In the arguments below, we shall anticipate over results to be fully 
demonstrated only a posteriori, via the formal developments in Sect. 4, and the explicit 
calculations in Sect. 5 and in Paper II (see also Refs. |10|, |1T], |1^, 0, p!9|]). 



The interactions between the various fields are described by the action + 5 A + a] 
To the order of interest, it is enough to preserve the linear coupling (5^4^(5 J" , where 



''^ ' - 5A^a(x) 



6^S 



6A^a{x)6A^,{y) 



(3.14) 

A+a 

a""'^^' 2 6A^a{x)SA'^{y)6A^,{z) 



A 



is the soft color current generated by the semi-fast quantum fluctuations. It is understood 
here that only the soft modes with ^ bA~^ are kept in the products of fields. In going 
from the first to the second line in this equation, we have expanded in powers of a^, and 
kept only terms which are linear or quadratic (Gaussian approximation). Note also that 
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we use compact notations, where repeated indices (coordinate variables) are understood 
to be summed (integrated) over. 

Given the separation of scales in the problem, we expect the color charge 5pa = SJ^ to 
be the "large component" of the current; that is, we expect the dominance of the eikonal 
coupling 6A~6pa (cf. Sect. 2.1). This will be indeed confirmed by the calculations, which 
show that the only current correlators to be retained to LLA are correlators of 6pa. 

Consider first the two-point function (SA^SA^). Quite generally, the following relation 
holds: 

{SA^^{x)SA^,{y)) = Jd'zJ d\GZci^,z) {6Jl{z)54{u))G^j:,,{u,y), (3.15) 

where Gr {Ga) is the retarded (advanced) propagator of the soft fields in the background 
of the tree- level fields and sources. To the order of interest, this can be computed in the 
mean-field approximation, that is, by inverting the following differential operator: 



Gj,{x,y)[A, p] - 



5A^'{x)5A''{y) 



(3.16) 



in the subspace of soft fields, in the LC gauge 5^4+ = 0, and with appropriate boundary 
conditions. Eq. ( p.l5| ) is simply the statement that the current- current correlator acts as 
a self-energy for the soft field two-point function (see Sect. 4.2 for a formal proof). To 
LLA, this equation can be simplified as follows: 

Among the various components of the polarization tensor {6J^6J'^), the logarithmic 
enhancement shows up only in the charge-charge correlator: 

Xabix,y) = {5paix)5pbiy)) . (3.17) 

Specifically, we shall see that x is a quantity of order ©(a^ ln(l/6)pp). 

To gain some more intuition, let us use as an example the contributions to Spa coming 
from the Yang-Mills piece of the action, Sym = I d'^x{—F^j^/4:) : 

6pa{x)\y^, = 2gr"^J^tiSK{x) + gr^'id^alixMix). (3.18) 

(The additional contributions to 6pa, coming from Sw, will be presented in Sect. 4.3.) 
The first term in the r.h.s. is linear in the quantum fluctuations a* and in the electric 
background field J-'~^^ (this is kinematically allowed because the semi-fast longitudinal 
momentum can flow from a* to JF+*). It is only this term that contributes to x to leading 
order in a<j. Indeed, it generates tree diagrams like the one depicted in Fig. |I|.a, where 
the internal line represents the propagator (a^a^) of the semi-fast gluons, and the soft 
external legs couple to the fields SA~. Physically, this diagram describes the emission of 
an on-shell, or real, semi-fast gluon by the classical source, which is a radiative correction 
to the direct emission of the soft field. This diagram will be evaluated in Sect. 5.2, but 
it is clear by power counting that it is of ©(a^pp). By contrast, one-loop diagrams which 



23 




a) b) 

Figure 1: Tree (a) and one-loop (b) contributions to x to quadratic order in p. Wavy 
lines with a blob denote insertions of the electric background field ; the thick internal 
wavy lines are propagators of the semi-fast gluons; the external dotted lines carry soft 
momenta, and couple to the fields 5A~ . 



are of the same order in p are necessarily of higher order in (7, and thus negligible for the 
present purposes. [An exemple of such a diagram, of 0{a1pp), is shown in Fig. |l].b; this 
is generated by the second piece, quadratic in a\ in the r.h.s. of eq. ( p.l8| ).] 



When computing Q{x,y), eq. ( p. 13 ), to LLA, the temporal non-locality in x can be 



ignored. To see this, consider eq. ( p.l5| ) with = and {5J^5J^) replaced by x. Since 
the external lines Gr and Ga are soft {k^ ^ &A+), the typical energy scale which controls 
the non-localities — z^ and x^—u^ (and therefore also —u^) is k~ = Q\/2k'^ ^ A^/6 
(cf. eq. ( p.llD ). This is much larger than the typical energy p~ ^ hr /h of the semi-fast 



gluons contributing to x- Therefore, one can neglect the evolution of 5p during the small 
time interval z'^ — ^ b/A~, and replace x{z~^ ~ '^~^) by its equal-time limit, which is 
independent of time: 

X{^,u) = xiz^ = u+,z,u) = a.ln^ x{z,u) . (3.19) 



As indicated above, it is only after taking the equal-time limit inside x that the logarithmic 
enhancement becomes manifest. The integrations over z~^ and then set p~ = in both 
external propagators, thus finally yielding 

{5A\x+,x)6A^{x+,y)) ^ f d^z f d^uG'-{x,z,p- = 0) x{z,u) G-^{u,y,p- = 0), (3.20) 



where we have omitted the subscripts R and A on the propagators since the boundary 
conditions in time become irrelevant as p~ = 0. In particular, since G = for p~ = 0, 
the two-point functions involving 6A^ (like {6A^6A^)) vanish to LLA. 

We now consider the induced mean field SA'^, which can be obtained by solving 
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eq. ( |3.5| ) to the order of interest: 



/ 6S 




\6Aii{x) 


A+SA+a) 



{a)=0 



6Aiiix)6Al{y) 



A 



(3.21) 



with an induced current: 



6^S 



^^"^^^ - 2 6A{i{x)6A'^{y)6A^{z) 



A 



<(y)a^(z)) + (6A^,{y)6A'Az 



(3.22) 



which involves two types of contributions: (a) A piece proportional to the 2-point function 
{a" a''') of the semi-fast gluons, and [b) a second piece involving the induced correlator 
{6A'^6A^) of the soft fields, as given by eq. ( |3.15| ). Let us explain this in more detail: 

(a) The first piece is the expectation value {SJ^) of the current (|3.14|) . Once again, 
only the /i = + component, i.e., the induced color charge density : 



[x) = {5pa{x)) 



(3.23) 



matters to LLA: a = ©(a^ ln(l/6)p). In Fig. ^.a we show a diagram contributing to a 
at lowest order in p; this is obtained by evaluating gf°''"^{{d~^al)al) to linear order in p 
(cf. eq. (|3.18|) ), and will be computed in Sect. 5.1. Note that the classical source p is 
represented here as a continuous line; this is to suggest the physical origin of p (namely, 
fast partons moving on straight line trajectories at the speed af light), and also the fact 
that its coupling to the semi-fast gluons is non-local in time (cf. Sect. 4.1). 

(6) To LLA, the second piece in eq. (|3.22| ) involves only the transverse field correlator 
{SA'SA^), which is proportional to x (cf- eq. ( |3.2CI| )). We denote this as: 



5Ju = -t: 



6^S 



2 dAt^dA^dAi 



A 



G'-xG 



-J 



(3.24) 



This is of order gx ~ gaslti{l/b)pp ; thus, in the saturation regime where gp ~ 1, it is as 
large as the induced charge a. 

By using Fig. |l].a for x, we obtain the contribution to 6J'^ depicted in Fig. |[b, which 
should be compared to Fig. ^.a for a : The logarithm ln(l/6) in the latter is generated 
by the integration over the internal p'^ momentum, which is restricted to the strip (|3.8|) . 
By contrast, in Fig. ^.b the loop momentum is soft, but the propagator {6A^6A^) is itself 
of order asln{l/b), since proportional to x- Thus, the logarithm is now generated by the 
subintegral giving x. Note that, while a describes the direct polarization of the semi-fast 
gluons, SJ''^ is rather an indirect effect of the latter, which induce a current after first 
modifying the propagator of the soft fields. 

To summarize, J'^ = 5^^(T + 6J''^ . As we shall verify in Sect. 3.3, where 6J'^ will be 
constructed explicitly, the current JT"^ is covariantly conserved, "DuJ^ = 0, as necessary for 



the mean field equation ( |3.21| ) to be well defined. (Here, = d^ — igA^ is the covariant 
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a) b) 

Figure 2: Contributions to the induced current to low orders in p. (a) A contribution to 
(7a ; the continuous line represents the classical source p. (b) A contribution to 5JI^ ; the 
thick internal line, with momentum p"*", is the propagator (a*a*) of the semi-fast gluons; 
the other, thiner, internal lines are the propagators and G^* of the soft fields. 



derivative constructed with the tree- level background field.) Since is independent of 
time, the LC-gauge solution to eq. ( |3.21|) is of the form (cf. Sect. 2.3): 



5At = 6A- = 0, 
SA'M = Jd'yG:i{x,y,p- = 0)J,\y). (3.25) 

At this stage, the gluon density in eq. ( |3.13D can be compactly written as follows: 

g{x,y)=A^,A^^+ {G^^J,),Al^ + A^,{G^^J,)g+ {G^-xG'%y. (3.26) 

The logarithmic enhancement is implicit in a (which enters J^) and In fact, to LLA, 
the induced charge (j{x) can be replaced in eq. ( |3.26| ) with a two-dimensional charge 
density localized at = : 

aa{x) 6{x')asln^(Ja{x^), (3.27) 



where 



as\n^aa{xi_) = J dx cra{x ,Xi_) (3.28) 



is the effective color charge density in the transverse plane. As shown by the last equality, 
the logarithmic enhancement becomes manifest only after integrating over x~. 

To understand this last operation, recall the relevant longitudinal scales: aa{x) is the 
color charge induced by the semi-fast gluons, and it spreads over longitudinal distances 
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1/A+ ^ ^ 1/6A+ (this will be verified by the explicit calculations in Sect. 5.1 and 
Paper II). Thus, cTq sits further away from the light-cone than the original source pa 
(which, we recall, has support at ^ a;~ ^ 1/ K^), but it is still very close to x~ = on 
the resolution scale of the soft fields 5 A''. In convolutions like 

{G^~a),^ j SyG'-{x,y)a{y), (3.29) 

the integral over y~ runs typically over large distances ^ 1/6A+, where it is indeed 
legitimate to replace a{x) as shown in eq. ( |3.27| ). 

Since eqs. ( |3.20|) and ( p.24| ) involve similar convolutions, there too one can replace 

Xab{x,y) -> S{x~)asln^Xab{x±,y_L)6{y'~), 
as In]- Xab{x±,y±) = f dx' [ dy' Xabix,y). (3.30) 



b 

In fact, X turns out to be even more localized than a : Indeed, x is explicitly proportional 
to the tree- level source p and/or electric field since it is generated by vertices which 
have this property (like the first term in the r.h.s. of eq. ( ^.18^ ; see also Fig. ^a). Thus, 



x{x,y) has support at ^ x ,y ^ 1/A+, like the original source p. 

The previous considerations show that the longitudinal and temporal structure of the 
quantum corrections do not matter when computing soft correlation functions: when 
"seen" by the soft fields, the (semi)fast gluons appear as a static color charge localized 
near the light-cone. This is consistent with the physical picture in Sect. 2.1, and also 
clarifies its limitations: This picture holds only to LLA, which is the accuracy to which we 
have a strong separation of scales. To this accuracy, only the one-point (a) and two-point 
(x) correlations of the quantum charge 6p must be kept (cf. Sect. 4.3 below). This is 
why, for matching purposes, it is enough to consider the two-point function (A^A^), as we 
do here. This is finally given by: 

g{x,y) ^ A^Al + asHl/b){{G'''J,),A}j + J^AG"'Ju)y + {G'-xG-')sy}, (3.31) 

where the logarithmic factor ln(l/6) is now explicit. After averaging over p with weight 
function Wa[p], eq. (|3.31|) gives the gluon density at the scale feA"*", as computed in Theory 
I to LLA. By construction, this must be the same as the gluon density in Theory II at 
tree-level. That is, the following equality must hold: 

(^t^;:),A = (^;:^;:)a + a.ln(l/6)((G^'^X)x^;^^ 
where, e.g., 

(^i4)A ^ Jvp Wa[p] A^M^aiy) , (3.33) 

together with a similar definition for {A^A'^)bA in terms of WbA- Eq- ( p.32D is the evolution 
equation for the gluon density, and can be used to derive the evolution iyA[p] — > W^bAM 
of the classical weight function, as it will be explained shortly. 
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To avoid cumbersome notations, in what follows we shall not always distinguish be- 
tween a and a, or between x and x\ these quantities differ just by their arguments and 
by a factor of asln(l/&) (cf. eqs. ( |3.27| )-( p.28| ) and ( |3.30|) ), which will be written down 
explicitly whenever needed. 

3.2.2 Tree-level calculation in Theory II 

In this subsection, the evolution equation for WV[p] will be established in two steps: 
First, we shall show that the quantum corrections described previously can be generated 
by adding a Gaussian "noise term" to the Yang Mills equation ( |2.2| ). Then, we shall 
absorb this noise term into a redefinition of the classical source and weight function. 

Our starting point is the following Yang Mills equation: 

[D„ = 6^+{pa{x) + Ua{x)), (3.34) 

which in addition to eq. ( p.2| ) contains a random source i^aix) (the "noise term") chosen 
such as to generate, via the solution to eq. ( p. 34 ), the same two-point function (AM*) as 



obtained in Theory I after including the quantum corrections (cf. eq. ( |3.31| )). Thus, Ua 
plays the role of the fluctuating charge 6pa of the semi-fast gluons. Based on this analogy, 
we take i>a to be static and have the same non-trivial correlators as 6pa, namely (with 
{■ ■ ■)u denoting the average over u) : 

{j^a{x))^ = (Taix), {Ua{x)Ub{y)) u = Xab{x,y). (3.35) 

From the previous subsection, we recall that the precise longitudinal structure of cr and x 
does not matter for the calculation of soft correlations. As we shall shortly discover, to a 
large extent this is also true for the evolution equation, which involves only the integrated 
densities aa{x±) and Xab{x±,y±), eqs. ( |3.28| ) and ( |3.3CI| ). Still, for this equation to be well 
defined, it will be important to recall that the induced charge aa{x) has a longitudinal 
support which is separated from that of the original source p (cf. the discussion after 
eq. ( p.28| )). This property is relevant here since the solution to the classical equation 
( ^.341 ) involves the color source via path-ordered exponentials in x~ , cf. Sect. 2.3. To 
keep trace of that, and for consistency with the first equation (|3.35| ), we choose I'aix) to 
have support in the range 1/A+ ^ x~ ^ l/feA"*", so like a. 

We first verify that the classical stochastic problem in eqs. ( p.34[ )-( |3.35D generates the 



same correlator {A^A^) as the quantum problem in Sect. 3.2.1. The solution y4'^[p, z/] to 
eq. ( p.34| ) can be read off Sect. 2.3: if A^lp] is the solution ( |2.18| ) to eq. (p. 21), then, clearly. 



Af'lp, u]{x) = S^^'A'Ip + z/](f). (3.36) 
For fixed p, this is a random variable whose correlations can be obtained by expanding 



in powers of u, and then averaging over u with the help of eq. ( |3.35| ). To the order of 
interest, we need the expansion of ( |3.36| ) to quadratic order (below, A^. = •A.l^[p]{x)) : 

1 5^Al 



6py 



Vv + 



^ 2 5py5p^ 



vyv, ^ A\\p\ + bA\\p,v\. (3.37) 
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Indeed, after averaging over z/, both the hnear and the quadratic term will contribute to 
order In (1/6), while the terms neglected in eq. (|3.37| ) would be of higher order. The 
classical two-point correlation function is then obtained as: 



where, to the order of interest, 

(ML). 



At^Ay + (ML).^; + ^L(K)- + ('^^LK)- (3-38) 



5py 



5pz 



1 



2 5py5pz 



Spu 



(3.39) 



These classical correlators coincide with the corresponding quantum corrections SA^, 
eq. ( |3.25| ), and {6A^6A^), eq. ( p.2q ), as we demonstrate now: 

For the two-point function, one can use the following identity, to be derived shortly. 



5^1 



5py 



[x,y,p 



0), 



(3.40) 



to conclude that the expression ( p.39| ) for {6A''6A^)i, is indeed the same as eq. (|3.2CI|) for 



{6A'6A^). Eq. ( ^.4U| ) can be proven by noticing that, to linear order in 6 A'' ^ u, 
which must be the same as the linear-order term in eq. ( |3.37| ). 

Eq. ( p.40|) also shows that the first term in the r.h.s. of eq. ( p.39| ) for {SA"^),^ is the same 
as G''~a, which is one of the contributions to 6A^ in eq. ( |3.25[ ). The other contribution, 
involving 5^7^^, can be similarly identified with the term proportional to x sq. (|3.39|) . 
More rapidly, the equality between {6A'')^ and 6A^ can be established by showing that 
they both satisfy the same equation of motion. In the case of {6A^)^, this is obtained by 
averaging eq. ( ^.34| ) over u, and can be rewritten as: 



'6S- 



YM 



6Aii 



A+SAj 



SfM-iPa + 0-a) 



(3.41) 



where Sym is the Yang-Mills action. For 6A\ the corresponding equation is rather ( |3.21| ), 
which is equivalent to: 



/ 6S 
\6M 



A+SAj 



S,,-cra 



(3.42) 



SinceQ A~ = 6A~ = 0, this is formally similar to eq. ( |3.41D . The final identification 
between these two equations then follows after recalling that the corresponding two-point 
functions {6A^6A^) (which enter both eq. ( |3.41 ) and eq. ( |3.42D ) are also the same. 

- like 



^■^By "SA = 0" , we mean here that {5 A ) = 0, and the two-point functions involving SA 
{SA^SA"^) — vanish as well (cf. the discussion after eq. ( 3.20 )). 
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To summarize, the classical two-point function ( p.38| ) coincide with the corresponding 
quantum function Q{x,y), eq. ( p.31|) . This allows us to write: 



(3.43) 



where the second equality is just a rewritting of the average over z/ as a functional integral 
with the Gaussian weight function: 



W[u;p] ^ e-^T^''^^exp|-i(z/-a).x;,J(^-^).} 



(3.44) 



(This depends upon p via x ci^nd a.) After also averaging over p, with weight function 
Wa[p], eq. (|3^) must be the same as the gluon density generated in Theory II at tree-level 
(cf. eq. ( |3.32|) ). This requires: 



■ Vp W,4p] Al[p]Al[p] = Jvp W4p] J Vv W[v ; p] A\\p + v\ A\\p + v\ , (3.45) 
which is satisfied provided 



WbA[p] 



VvWaIp - v]W{v \ p - v] . 



(3.46) 



This is a functional recurrency formula for W^[p], whose r.h.s. can be expanded as follows: 

WA[p-v]mv^p-v] - 



1 - + 7:y.yy^^]wA[p\W[v,p\, (3.47) 

bp^ 2 bpjpy 



and then integrated over z/ to give, to the order of interest. 



(3.48) 



A priori, the convolutions in the r.h.s. of eq. ( p. 48 ) involve three-dimensional integrals, 
like 



Opx 



6 



5pa{x) 



WAOaix) 



(3.49) 



Recall, however, that a is distributed over distances 1/A+ ^ ^ l/feA"*", and that the 
logarithmic enhancement is seen only after integration over x~, cf. eq. ( p.28| ). In the limit 
6 — > 1, in which we are eventually interested, eq. ( p.49|) can be replaced with 



7^ Wkf^x] = cts In 7 / d^x^ - — — ^ 

opx J dpa{x^,x±) 



WAaa{x±) 



(3.50) 



where the functional derivative is to be evaluated at = 1/A^. 
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It is convenient to introduce the rapidity variable^ r = ln(P+/A+) = ln(l/x), where 
is the total longitudinal momentum of the hadron, and x is Bjorken's x (cf. Sect. 
2.1). Then, ln(P+/6A+) = r + Ar, with Ar = ln(l/fe). After also renoting Wa = Wr, 
WbA = Wr+Ar, and a;^ = 1/A+ = x~ , eq. ( |3.48|) is rewritten as 

W.,.M - WM = .,Ar {1^_^J!_^|„.;,,,| - ^|H->j} . (3.51) 

where Pt{x±) = p{x~,x±), and the convolutions are now to be understood as two dimen- 
sional integrals, cf. eq. (|3.50|) . For instance: 

[WrCT^] = I (fx^——^ -\Wraa{x^)\. (3.52) 



5pr{x) J 6pa{x^,X± 

According to eqs. ( p.51|) -( p.52| ), the evolution from Wr[p] to Wt+AtIp] is due to changes 
in p within the rapidity interval (r, r + Ar), where the quantum corrections are located. 

Note that the variable r in the above equations acts simultaneously as a momentum 
rapidity [r = ln(P"'"/A"'")], and as a space-time rapidity [r = \ia{x~ /xq), with Xq some 
arbitrary scale of reference; e.g., Xq = 1/P^]. This is so since the support of the quantum 
corrections is correlated to the longitudinal momenta of the modes which have been 
integrated out. These equations suggest a space-time picture of the quantum evolution 
where the color source is built up by adding contributions in successive layers of (space- 
time) rapidity, from x~ = up. 

By taking the limit Ar = ln(l/6) — > 0, we finally obtain a functional differential 
equation for the evolution of the weight function with r = ln(l/a;): 



This renormalization group equation (RGE) has been originally derived in Ref. |TT 
although its longitudinal structure — i.e., the fact that the functional differentiations are 
to be taken with respect to the color source in the highest bin of rapidity (r, r + dr) — 
has not been recognized there. (In Ref. [|ll] the differentiations are rather taken w.r.t. 
the integrated charge density p{x±) = J dx~ p{x).) To understand the importance of this 
structure, recall from Sect. 2.3 that the classical solution A^[p] — and therefore also the 
quantum corrections a and x? which are functionals of p2| — depends upon p via 



Wilson lines in x (cf. eq. ( |2.29|) ). Thus, in evaluating eq. (|3.53|) , we shall be led to 



differentiate 132 



V\y^) = Pexp|«(7 f^^^dz- j d^z^{y^\^\z^) pa{z- ,z^)tA , (3.54) 



^^We prefer to denote this as r, rather than the more conventional y, to avoid possible confusion with 
the coordinate variable y, which also enters the subsequent equations. 
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with respect to Pt{x). Note the upper hmit x~ in the integral above: this occurs since 
the original source p has support up to x~ ~ x~. Because of the path-ordering in x~ , the 
functional derivative of eq. (|3.54|) is as simple as a normal derivative: 

■V^iy^) = tg{x^\-^\y^)T'^V\y^). (3.55) 



6pa{x^,X±] 

By contrast, its derivative w.r.t. the integrated charge Pa{x±) is not well defined. 

Note also that the color source which appears in eqs. ( |3.54D -( ^.55| ) is the covariant- 
gauge source pa, and not the LC-gauge source pa- It is only with respect to pa that we 
know the classical solution explicitly. This brings us to the necessity to rewrite the RGE 
( |3.53| ) as an evolution equation for WV[p], which will shall do in the next subsection. 



Let us conclude this subsection with a few remarks on eq. ( p.53|) . This can be rec- 



ognized as a functional Fokker-Planck equation |]3^, with r playing the role of time. 
It describes diffusion in the functional space spanned by p, with (p-dependent) "drift 
velocity" agcr and "diffusion constant" asX- (Iii this language, the recurrency formula 



( |3.46|) can be viewed as a functional Chapman- Kolmogorov equation [0.) Alternatively 



eq. ( p.53| ) is like a functional Schrodinger equation in imaginary "time" r. It can be trans- 



formed into an hierarchy of ordinary differential equations for the correlation functions 
{pp- ■ ■ p)r generated by Wr[p] ■ For instance, by multiplying it with p{x)p{y) and func- 
tionally integrating over p, one obtains an evolution equation for the two-point function: 

^ {paix)pbiy))^ = as (S{x~ - x~)aaix:i)pb{y) + 5{y~ - x-)pa{x)ab{y_i) 

+ 6{x~ -x:;)6{y' -x~)xah{x^,yi_))^, (3.56) 

where (■ ■ ■)r denotes the average over p with weight function W^r[p] • For strong fields and 
sources, that is, in the saturation regime, this equation involves also the higher n-point 
correlators (n > 3), since a and x ^^re non-linear in p But in the weak source limit. 



where the non-hnear effects can be neglected, eq. ( p.56|) closes at the level of the two-point 
functions, and becomes equivalent to the BFKL equation, as demonstrated in Ref. [|l3 
(see also Sect. 5 below, and Paper II, where the BFKL equation will be recovered as the 
weak field limit of the general evolution equations). 



3.3 Evolution equation with a covariant-gauge source 

For the RGE ( |3.53| ) to be of practical use, its coefficients a and x must be known exphcitly 
in terms of the tree-level source p. This requires, in particular, an explicit expression for 
the classical solution A^[p\ (since a and x ^-^e functionals of A^] cf. Sects. 4, 5 and Paper 
II). From Sect. 2.3, we know that such an explicit expression can be obtained if the light- 
cone gauge solution is expressed in terms of the covariant gauge (COV-gauge) charge 
density p. This observation, together with the gauge invariance of the weight function 
(W/a[p] = Wa[p]), makes it more convenient to use p as the classical source variable to be 
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averaged over. In this subsection, we shall rewrite the RGE in this new variable. This is 
non-trivial since, as it will be explained shortly, the gauge rotation from p to p is itself 
subjected to quantum evolution. 

Recall indeed that the definition of the COV-gauge depends upon the classical color 
source in the problem: For a given static source pa{x) in the LC-gauge, the COV-gauge is 
the gauge where the classical field generated by p has just a plus component, = 6'^~^a, 
with a hnearly related to p (the COV-gauge version of p). Thus, if the classical source 
is modified by quantum corrections (say, from p to p + in the LC gauge), then what 
we mean by "the covariant gauge" changes as well: This is now the gauge where the new 
classical field, as generated by the total charge p + i', has just a plus component, say 

_ ^M+Q,^ with a linearly related to the COV-gauge source p + u. [Note that we denote 
with a tilde (a bar) quantities in the COV-gauge associated with p (respectively, p + u).] 
Clearly, the rotation between the old and the new COV-gauges depends on z/, and thus 
upon the quantum corrections. In what follows, we shall construct this rotation explicitly. 

Recall first the classical solution corresponding to p, as constructed in Sect. 2.3. This 
is summarized in the following equations: 

A'[p]{x) = -^U{x)d'U\x), 

U\x~,Xi_) = Pexp|i(7 J dz~ a{z~ ,x±)^ , 
-Vla(f) = p(f), (3.57) 

where W{x) is the gauge rotation from the LC-gauge to the original COV-gauge (cf. 
eq. m^). 

When the classical source changes from p to p + z/, the corresponding solution is 
obtained by replacing p p + z/ in eqs. ( p.57| ) : 

A'[p + iy]ix) = -U(x)d'uUx), 
9 

W{x^,x±) = Pexpji^f J dz~ a{z~ , x±)^ , 

-Via = p + u. (3.58) 

We write W = 6W t/^, with 6U'^ the gauge rotation from the original COV-gauge towards 
the new one. It turns out that this has a relatively simple expression in terms of u. 



To see this, we set a = a + 6a, with a defined by the last equation ( p.57| ) and Sa 



proportional to u. These two fields have no overlap in the longitudinal direction: a, 
like p, has support at ^ x~ ^ while z/, and therefore 6a, is rather localized at 

1/A+ ^ x~ ^ 1/6A''". Then, the path-ordered exponential in eq. (p.58|) simply factorizes: 

W{x',x^) = 6U^{x)U^{x) ^ exp[ige{x- - x;)6a{x^)}u^{x) , (3.59) 
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where = l/A"*", Sa{x±) = J dx Sa{x), and the approximate equahty holds as 6 — > 1. 
In particular, within the support of p, W reduces to [/^, and therefore 



P 



U^pU = U^pU = p. (3.60) 



Thus, remarkably, the original source p is not affected by the supplementary gauge rotation 
6W induced by the quantum corrections. By using this property, and subtracting the 
equation satisfied by a from the equation for d, we deduce 

- \/l5a{x) = u{x) , (3.61) 

which shows that Sa is linear in u. 

We are now in a position to reformulate the RGE by using the COV-gauge source as 
a variable. The analog of eq. (p.45|) reads: 

I Vp W,a[p] A'M^llp] = jvp W^[p\ I Vu W[u ; p] ^^[p + u] ^;[p + u] . (3.62) 

Note the assignments of tilde's and bar's in the r.h.s.: This is what we need since, first, 
p is the COV-gauge source at tree-level, and, second, the classical solution corresponding 
to p + z/ is the same functional of p + z/ as the tree-level field in terms of p (cf. eqs. ( |3.57| ) 
and ( |3.58| )). By also using p = p (cf. eq. ( |3.60| )), the r.h.s. can be rewritten as: 

JvpWa[p] JvuW[u;p]A',[p + u]Ai[p + u], (3.63) 

which leads to the following recurrency formula (compare to eq. ( |3.46|) ): 

WtA[p] = JvuWaIp- ^] >V[z> ; p - z7] . (3.64) 

Thus, in this framework, it is natural to define the noise correlations directly for u. We 
write (compare to eq. ( PTH 



(z7"(f)), = a:,(f), {u'^ixyiy)), = rj{x,y), (3.65) 

which together with eq. (|3.64|) leads to a RGE for H^r[p] with the same formal structure 
as eq. p.53| ), but with p ^ p, X Xci, and ad (Xci. It just remains to specify the 
functions aci and Xd in eq. ( p.65| ). 

These can be obtained from the corresponding correlations in the LC gauge, eq. ( |3.35| ), 
by an appropriate gauge rotation. We have: 

V = tJ^vU = 6U^[WuU)6U = 5U^v5U , (3.66) 

with the correlations of z/ following trivially from eq. (|3.35| ): 

{i^a{x))u = aa{x) = Ulf^{x) ab{x), 
{ua{.x)h{.y))u = Xabix,y) = Ul^{x) Xcd{.x,y) Udh^v) ■ (3.67) 
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But the relation between v and v is non-linear, since the "small" gauge rotation bV^ is 



itself dependent on v. To the order of interest, one can expand eq. (|3.66|) to quadratic 
order in z/, to obtain (cf. eq. ( p.59 )) : 



6{x — x^)6a{x_i_), V = u + g -[6a{x±), u], 



(3.68) 



where 6a is to be taken to linear order in u, and the factor 1/2 occurs since, as 6 — > 1, 
z>(x) ^ 6{x~ — x~)h'{x±) and 9{x)6{x) = (l/2)5(x). By using the equations above, and 
after simple algebra, one finally obtains 



a 



d 

9 nabc 



6a''^{x) ^ / d'y X6c(f , y) {vA ^ . 



(3.69) 
(3.70) 



The correction —5aci is a "counterterm" which takes care of the spurious classical polar- 
ization which appears in the covariant gauge (cf. eq. ( |3.77D below). Note that 5a ci ~ 
5fa;jln(l/6)pp, which in the saturation regime {gp ~ 1) is as large as a. 

To summarize, in addition to the straightforward gauge rotations ( 3.67| ), the RGE for 
Wr[p\ differs from the corresponding equation ( p.53| ) for Wr[p\ also by the correction ( |3.70| ) 
to the induced charge density. This latter is the consequence of the quantum evolution of 
the "covariant gauge" itself. 

By exploiting the previous results, it is now straightforward to derive explicit expres- 
sions for the classical correlations generated by z/ (cf. eq. (^.39|)), to be then compared 



with the quantum correlations in Sect. 3.2.1. This will allow us to understand the struc- 
ture of these correlations in more detail. To this aim, one needs the expansion of ^*[p + z/] 
to quadratic order in z/, which reads (cf. eqs. ( |3.58|) and ( p.59| )): 



5A' 



-U fe-^3^" 9^e^^^") f/^ = U6A'U\ 
9 ^ ' 



&'ba + g 



1 r 



5a, d^ba 



(3.71) 
(3.72) 



Here, ba = Q{x~ — x~)(5a(x^), so that the induced field bA" has support only at positive 
and (relatively) large x~ <^ l/A"*" (compare with the original field in eq. ( p.31| )). 
This property is the consequence of the fact that the induced source itself has support at 
x~ ^ l/A"*", together with the retarded boundary conditions used when solving eq. ( |3.34|) . 



Note the relatively simple expression for the linear term in eq. ( p.72| ) for SA^. This is 
the same as (cf. eq. ( |3.61| )) 



6A\x) 



linear 



9{x- -x-)[^v 



d'yGl{x,y,p- = 0)u{y), (3.73) 
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where Gq (P^ = 0) = —p^ip^PjJ, with l/p'^ = + ie) , is the free LC-gauge prop- 

agator, cf. eqs. (|2.8|) and (|3.84|) . A priori, the emergence of the free propagator might 
seem surprising; one would rather expect the following Green's function (with p~ = and 
G'*"[^*] as introduced in eq. (|3.15| )) : 



(3.74) 



which is the propagator of the LC-gauge fluctuations 6A^ in the presence of the COV- 
gauge background field^ = In fact, eqs. (|3.73|) and ( p. 741) are consistent each 

other: indeed, a{x) is localized near x~ = 0, while the integral over y~ in eq. ( 3.73 ) runs 
from y~ = x~ > to x~ > y~ . That is, the whole propagation takes place far outside the 
support of a, where G*^ = G}^ indeed (cf. Sect. 6). 

The correlations {SA'^)^, and {6A'^6A'^)iy are simply gauge rotations of the corresponding 
correlations of SA^ (cf. eq. (|3.71|) ). To the order of interest, one immediately obtains: 



SA:{x)SAI 



(fz / (fu{x_L\ 



vi 



\yi.) 



(3.75) 



which, as expected, is the (gauge rotation of the) quantum correlator ( |3.2CI|) . 
The one-point function {6A^)iy is more interesting. It reads (for x" > x^) : 

(9* _ i 
2 



6a, d^6a 



(3.76) 



where the second term in the r.h.s. is proportional to x- After simple algebra, one finds 



6a, d^6a 



d' 1 



6T 



(3.77) 



where 6aci is defined in eq. ( p.70| ), and the vector current 



6Ji = -r'\d'd^ -S'^Vl){6a,&6a,), 



f-bc^(^^QiQj _ ^ij^2^^^^^^ d^6ac)^ (3.78) 



is transverse, 9^6 Jl^ = 0, as required by gauge symmetry. (In the background field LC 
gauge, this becomes 'D'^6J'^ = 0, with "D* = 9* — igA\) By the matching conditions, this 
is the same as the current ( p.24|) generated by soft fluctuations in Theory I. In particular, 
we note that 6J'^ = : no charge density is induced by these fiuctuationsQ. 

The induced mean field reads finally (cf. eqs. ( p.76| ), ( p.77| ) and (|3.69| )) 



{6A% = ^a + ^6J^ 



(3.79) 



^^Note that the fluctu ations transform homogeneously under the gauge rotatfons of the background 
field; see, e.g., eq. ( 3.71 ). 
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This could be also verified via a direct calculation of dj+ = gf'''"'{{d+5Al)SAl,) (cf. eq. ( [4.421) ) 
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and coincides, as expected, with the (gauge rotation of the) quantum correction SA\ 
eq. (|3:25|) . 

Note that, unhke the tree-level field A^, the total mean field including quantum correc- 
tions, + 6A\ is not a "pure gauge" anylonger (because the current JT,y has non-trivial 
transverse components J'i, cf. eq. ( p.78| )). Rather, it is the fluctuating field + i^], i.e. 
the solution to the classical equations ( |3.34]) for a fixed value of z/, which is a pure gauge 
(cf. eq. (|3.58| )). This is as expected, since, in the classical effective theory defined by p, 
it is the solution A^[p] at fixed p which has this property. 



3.4 RG evolution in p 

Previously, the "semi-fast gluons" (i.e., the quantum fiuctuations to be integrated over in 
one step of the renormalization group analysis) have been defined as the fiuctuations with 
longitudinal momenta inside the strip < \p'^\ < A"*". But we have also mentioned that 
the leading logarithmic effects involve nearly on-shell quanta {2p^p~ k, p^), for which the 
strip restriction on implies a similar restriction on p~: 

hr < \p-\ < A-/b, (3.80) 

with A" = QI/2A+, and generic Q± (cf. eq. ( p.ll|) ). Thus, in practical calculations, the 
strip restriction can be imposed either on p"*", or on p~ . 

In Sect. 5, we shall perform calculations in the weak source limit and show that, 
to LLA, both strip restrictions lead indeed to the same results (namely, to the BFKL 
equation). For general, non-linear, calculations however, it turns out that it is more 
convenient to work with a strip restriction on p~ . The reason is that the tree- level fields 
and sources are independent of time, so any constraint on the energy of the quantum fields 
will be automatically preserved during their propagation. By contrast, the inhomogeneity 
of the background fields in x~ makes it difficult to maintain a similar restriction on p"*". 

To be more specific, let us note that the key ingredient in the computation of cr and 
X is the 2-point function {a^{x)a"{y)) of the semi-fast gluons in the background of the 
tree-level fields and sources A^ and p. For instance, eqs. (|3.14|) and ( p.23| ) imply: 



.... 1 

Gaix) = 



2 6A-ix)6Atiy)6A^M 



KivKiz)) . (3.81) 

A 



As we shall see in Sect. 4, this 2-point function can be obtained from the following 
time-ordered (or Feynman) propagator: 

iGZ{x,y)[A, p] = (T<(x)a;;(y)) 

= Z^l^^ (\a5{a+) <(x)<(i/)e^^°[^''^l , (3.82) 
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where, as compared to eq. ( p.2| ), the functional integral is now restricted to semi-fast 
fluctuations, and the action is evaluated in the Gaussian approximation: 



S[A + a, p] 



(fx / (fy a'^{x) 



5A^a{x)5Al{y) 



A 



O'biy) = So[A,a] 



(3.83) 



[The first two terms in this expansion vanish since A^ is a solution of the classical equations 
of motion SS/SA'^ = 0, and, moreover, the action itself vanishes on this solution: S[A, p] = 
; cf. Sect. 4.1.] Thus, this propagator is obtained by inverting the differential operator 
( p.l6|) in the LC gauge, and in the subspace of fields with momenta restricted to the strip. 

For sufficiently weak fields and sources, the propagator can be computed via an expan- 
sion in powers of A^ and p. In particular, when p ^ 0, and thus ^* ^ 0, it must reduce 
to the free propagator given in eq. ( 3.84] ) below. But in the saturation regime, where the 
background fields are as strong as p ^ 1/g and A"^ ~ 1/g, this propagator must be known 
exactly (i.e., to all orders in A^ and p). In that case, it is easier to work with a strip 
restriction on p~ : then, the p'^ momenta are unrestricted, which makes it convenient to 
use the x~ -representation and thus exploit the specific geometry of the problem (namely, 
the fact that the source p is localized near x~ = 0) in order to construct the propagator. 
This construction will be presented in Sect. 6. 

To conclude, from now on we shall impose the strip restriction on p~ , cf. eq. ( |3.80| ). 
This is consistent with the background-field gauge symmetry, since all the gauge transfor- 
mations of the background fields to be performed (like the rotations between the LC-gauge 
and the COV-gauge) are time-independent, and thus cannot change p^ . But this leaves 
the place for a potential ambiguity in the quantum gauge condition, as associated with 
the axial pole at p"^ = in the LC-gauge propagator. That is, when inverting eq. (|3.16|) , 
one has to specify not only the ie prescription at the physical pole for on-shell excita- 
tions (here, the Feynman prescription), but also the prescription for the unphysical pole 
in l/p"*". The latter is, however, fixed by our choice of the boundary conditions for the 
classical solution in Sect. 2.3, cf. eq. ( ^.57| ): For this to be consistent with the quantum 
evolution, one must use the "retarded" prescription l/{p'^ + ie) in the G*~ component 
of the propagator (see the discussion after eq. ( |6.33D ). The prescriptions in the other 
components then follow by hermiticity. 

In particular, the free propagator reads as follows: 



Gl\p) 



p 



p^ + le 
6'^Goip), 



Goip), G-,\p) 



p 



p^ — le 
2p- 



-Go{p), 



G^-{p)= PV^Go(p) 



pi 



(3.84) 



where Gq{p) = l/{2p^p — p]_ + ie), and PV denotes the principal value prescription: 

1 1 \ 



+ 



2 \ p+ — ie p+ + ie 



(3.85) 
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Other prescriptions to be later referred to are the "advanced" prescription previously 
used by Mueller and Qiu , and by Kovchegov and Mueller H — this is obtained by 
changing the sign of ie for the axial poles in eq. ( p.84 ) — , and the PV-prescription, for 
which l/p"^ = PV{l/p^) in all the components of the propagator. 



4 Quantum dynamics on a complex-time contour 

In this section, we describe the action S[A, p] which enters the quantum McLerran- 
Venugopalan model, eq. (|3.1| ), and derive Feynman rules for the computation of cr and 
X- Because of the temporal non-locahty of the classical equations of motion ( p.l7|) , this 



action, and the associated quantum dynamics, must be formulated along a contour in the 
complex-time plane. It turns out, however, that in the approximations of interest the 
contour structure is not essential, and one can restrict oneself to the dynamics in real 



time, as done a priori in Refs. [TO, ITT . 



4.1 The contour action 



By definition, the action S[A, p] must be such as to reproduce the classical equations of 
motion ( |2.17| ) at the tree-level. Thus, S[A, p] has the following structure: 



S[A,p] = - J d'xh^^Fr + Sw[A-,p], 
with 5*14/ [^^,p] constrained by (cf. eqs. ( ^.15| ) and ( |2.16| )) 

-J+(x) = -^Tt{p{x)W\x)T^W{x)} 



6S^ 



6A~ix) 



(4.1) 



(4.2) 



However, with the retarded current J~^{x) in eqs. ( |2.15| )-( ^.16| ), the constraint ( [4.2| ) can- 
not be satisfied, for any action. For instance, eq. ( [4.2| ) implies the following symmetry 
property, or Maxwell relation : 



5J+(x) 



w 



6A^iy) 6A-ix)6A^iy) M-(x)' 
which is however inconsistent with the current (|2.15| ), which rather yields 



(4.3) 



SJ^ix) 



^e(x 



y+)6^'\x-y)TT{p{x)[T^,T']} 
-gr'^p^[x)e{x+ -y+)5^'\x-y) ^ nf,(x,y). 



(4.4) 



That is, the two-point polarization function 11^, and all the higher vertices generated by 
J^{x), are necessarily retarded with respect to a;"*", while the corresponding amplitudes 
generated by Sw should be rather symmetrical (e.g., Iiab{x,y) = Ilba{y,x)). 
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Figure 3: Complex-time contour for the quantum MV model: C = C+ U C_. 



Thus, the non-local equation ( p.l7| ) cannot be generated from an action. This reflects 



the fact that this is an effective equation, obtained after some degrees of freedom (the fast 
partons) have been integrated out, with retarded boundary conditions, to generate the 
non-local current in the r.h.s. There is a standard procedure to circumvent this difficulty 
and construct a quantum version of the non-local theory (see, e.g., Ref. |^0[ and Refs. 
therein): This involves a generalization of eq. (|2.17|) to complex time variables which take 
values along the Schwinger-Keldysh contour depicted in Fig. ^. 

We call z the (complex) time variable along the contour, and reserve the notation 
for real times. The contour C may be seen as the juxtaposition of two pieces: C = C+UC_. 
On C+, z = takes all the real values between Xq and x~f (eventually, we let Xq — — oo 
and Xj oo). On C_, z = x~^ — it], with rj 0^ and x^ running backward from Xj 
to Xq. We define a contour 6'-function 9c : Oc{zi,Z2) = 1 if 2:1 is further than Z2 along 
the contour (we write then zi >- Z2), while Oc{zi,Z2) = if the opposite situation holds 
{zi -< Z2). We can formalize this by introducing a real parameter u which is continously 
increasing along the contour; then, the contour C is specified by a function z{u), and 
Oc{zi, Z2) = 6{ui — U2). The contour 5-function is defined by: 




5c{zi, Z2) = \ ^ ] 5{ui-U2). (4.5) 



The quantum extension of the MV model is a quantum field theory living on this contour. 
Specifically, we allow the quantum fiuctuations of the fields to take arbitrary values 
on both sides of the contour, and define a contour "evolution operator" by (for Z2>~ Zi) 



x) = Tc exp I i(7 J dzA (z,x)^, (4.6) 



40 



where J^^ is the integral running along C from Zi to Z2, and the operator Tc orders the 
color matrices A~{z) from right to left in increasing sequence of their u arguments. (Note 
that the ordering along the contour coincides with the chronological ordering on C+, and 
with antichronological ordering on C_.) In particular, for Zi = Xq and Z2 = Xq — irj (the 
end points of the contour), we write: 

Wc{x) = Tc expf^ig J^dzA-{z,x)Y (4.7) 

with running all the way along the contour. 

We are now in a position to present the quantum action for the non-local effective 
theory (below, d^x = Jcdz J d^x) : 

S[A,p\ = - J^d'x^F;,Fr + j d?Si:i[p{x)Wc[A~]{x)} = Sym + Sw (4.8) 

This is real, as it should, since the field A~ in the Wilson line ([4.7]) is in the adjoint 



representation: A = A^T'^. In the saddle point approximation, eq. ([4.8|) generates 



eq. (|2.17|) , as we verify now. Note first that (with Xq —00) : 

J^^^ = -j^^Pim-oo-^r,A^)T^W.,.^ix)}. (4.9) 

The tree-level field A'^[z) (the solution to 5S/SA^{z) = 0) takes identical values on both 
sides of the contour: A'^{x~^) = A'^{x~^ — irj) for any real x~^. For such a field, eq. (^^ ) 
reduces to the current J~^{x) in eq. (|4.2| ), for z irrespectively on C+ or on C_. To verify 
this, consider, e.g., z = x^ E C+; then, Wx+-oo = W{x'^) (cf. eq. (|2.16|) ), and 



^^—00— iri,x+ ^^—00— iri,x+ —irj^^x+ —iri,oo— iri^^(X),x+ : i.'^AO^ 



where we have used the group property of the contour operators ( [4.6|) . Since the tree- 
level field A~ is the same on both sides of the contour, the forward evolution in time from 
x~^ to 00 is compensated by the backward evolution in time from 00 — irj to x~^ — it] : 
M^x+-i77,oo-ir? W^oo,x+ = 1- Similarly (with T denoting anti-chronological ordering). 



^^—oo—iri,x+—iri 



W^{x+) = Texpl-ig£^ dz+ A- {z+ , x)\ . (4.11) 



Thus, at tree-level, and for z = x+ e C+, W.^.i^^zT" W,^_^ W^{x+)T''W{x+), in 
agreement with eq. ( [4.2| ). A similar conclusion holds for z = x^ — irj & C_. In both cases, 
the saddle point equation 5SI5A^^{z) = turns out to be equivalent to eq. ( p.lTp . 

The non-local piece of the action Sw describes the propagation of the color source p 
(= the fast degrees of freedom of the hadron) at the speed of light, in the background of 
the soft fields A^ (the classical field radiated by p plus the soft quantum fluctuations), 
and in the eikonal approximation. This gives rise to vertices of the type p(yl~)", n > 1, 
which are non-local in time (see Fig. ^). 
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Figure 4: A typical n-point vertex (here, n = 3) generated by the action Sw- This is 
non-local in time, but local in the spatial coordinates. The continuous line represents the 
source p, while the wavy lines are A~ gluon fields. 



Sy/ is invariant under the periodic gauge transformations VL{x) G SU(A^) which satisfy 

fi(-oo-i?7,f) = fi(-oo,f) (4.12) 
for any x. Indeed, under such a transformation, 

A'^(x) — ^ Vl{x)[A^' + {i/g)d^')n\x) , 

p{x) — > Vl{—oo, x) p{x)VL^ {—oo,x) , 
Wc{x) — ^ n{-oo-i7],x)Wc{x)n\ —oo,x), (4.13) 



and Sw is invariant because of the property ([4.12|) . Here, we are mostly interested in time- 
independent gauge transformations, which preserve the static character of the background 
fields, and for which the periodicity condition (|4.12|) is trivially satisfied. 

The transformation law for p{x) in eq. (|4.13| ) refiects its interpretation as the initial 
value of the color charge at Xq = —oo (cf. Sect. 2.3)). In fact, Sw can be also written 
as follows|3 : 

Sw = -j^ J^d^xTT{A-{z,x)W,,^^{x)p{x)}, (4.14) 

which emphasizes the fact that this is a gauge-invariant generalization of the linear vertex 
Jd^xpaA-. 



4.2 Contour Green's functions 

The generating functional for contour Green's functions is defined as in eq. ( p.2|) where 
the time variables now run along the contour, and the external source is extended to 



17 



To see this, use the identity idzWz-oc = gA {z)Wz~oo to perform the integral over z in eq. (|4.14| ) 
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a function jl^{z,x) on C. For instance, the (connected) contour two-point function is 
obtained as (with color indices omitted): 



^jf^{zi)Sj^{z2) 



j=0 



ec{Zl,Z2)G>^"'{z,,Z2) + ec{z2,z{)G<^''{Zi,Z2), 



(4.15) 



where, to hghten the notations, we have indicated only the time variables. In the second 
line above, we have introduced the two-point correlation, or Wightman, functions: 



>r{z,,z2)^{A>:{zMi{z2))c, 



■:r{^i,^2)^{Al{z2)A^Az,)),, 



(4.16) 



which, unlike the contour propagator Gq" {zi, Z2), are continuous at zi = Z2. These are 
the functions which enter the quantum corrections a and X] for instance, the contribution 
of eq. ( p.l8| ) to a reads (see also Sect. 4.3 below): 



cra{x)\ 



YM 



gr'%{al{y)al{x))\^^^ = gr'%Gr{x,y)\ 



y=x 



(4.17) 



For real time variables, the functions G^ and G^ can be used to construct the retarded 
(Gr) and advanced (Ga) propagators, which will also be needed: 

GR{x,y) = z9{x+-y+)[G>{x,y) - G<{x,y)], 

GA{x,y) = -te{y+ -x^)[G>{x,y) - G<{x,y)]. (4.18) 

As a simple example, let us construct the contour propagator for a free scalar field 
(the generalization to a free gluon field in the LC gauge is straightforward). For both time 
arguments x~^ and y'^ on C+, this coincides with the usual Feynman propagator, which 
reads (with p = (p+, p±)) : 



iGo{x^ - y^,p) 



dp 



-ip J/+) 



2tt 2p+p^ — p]_ + ie 

= \e{p+)e{x+ -y+)-e{-p+)e{y+ -x+)]e~"^ 

2p+ ^ ^ 

from which one can identify the free Wightman functions: 



(4.19) 



ix+-y+) 



G>ix+,p) 



■e 2p-t 



Q(x+,p) 



2p+ ' 2p^ 

or, in momentum space (with = 2p'^p~ — p\), 

G>{p) = 2ne{p+)6{p'), G<{p) = 27Te{-p+)6{p') 



e 2p-t 



(4.20) 



(4.21) 



The free contour propagator Gq is finally obtained by inserting Gq and Gq in eq. ( |4.15| ). 
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Returning to the gauge theory with action ( ^.Sp , the corresponding contour propagator 
Gc satisfies the Dyson-Schwinger equation with time arguments on C: 



5^S 



c 



5Ai'{x)5A''{z) 



A 



(4.22) 



Here, 6Il^j^{x,y) is the (contour) self-energy, and describes quantum corrections, and 



6A^'{x)6A•'{y) 
where we have denoted 



1^fiu{x)6c{x,y) + 6^-6u-Il'^{x,y) 



A 



nS,(a;,y) 



^IV 1 



w 



6A-{x)6A, (y) 



A 



(4.23) 



(4.24) 



In these equations, = — igA'^, and A'^ is the background field, to be eventually 
identified with the classical field in Sect. 2.3 (hence the notation). In particular, A^ = 0, 
and a second differentiation in eq. ( [4. 91) yields (with x^ , y~^ ^ C): 



I r'^Pcix) 6^'\x - y){9cix+, y+) - e^y^, x+)} 



= 0c(x+,y+)n>(x,i/) + ec{y^,x+)Yi<,{x,y), 
with time-independent 11^ and 11^ : 

n>,(x,i/) = -{g/2)r''^p,{x)5^^\x-y) = -U<ix,y). 
The contour self-energy admits a similar decomposition: 

The Dyson-Schwinger equation is then conveniently rewritten as 

V,,ix)G''c\x,y) + / d'zJ:^(x,z)G^c\^,y) = Si5c{^,y), 



c 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



with the total self-energy E 



c 

flU 



-'fi— 



By choosing x~^ G C+ and y~^ G C_ in eq. ( [4.28[ ), and using the decompositions ( |4.15| ) 
and (|4.27|) , we obtain, after simple manipulations, an equation for G^(x, y) in real time : 



V/{x)G<,{x,y) = Jd'z{i:nG< + ^<GA),Ax,y). 



(4.29) 



[The retarded (S/j) and advanced (S^) self-energies are defined in terms of the Wightman 
self-energies and as in eq. ( [4.18| ).] One similarly derives an equation for G-^{x,y), 
as well as the following equation for Gji{x,y): 



V%{x)G'^{x,y) - / d'ziJ:nGnrix,y) = g^^5^'\x-y). 



(4.30) 
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The two equations above imply a relation between and : 

G^ui^.y) = J d^zd\{GR{x,z)J:<{z,u)GA{u,y))^^. (4.31) 
A similar relation holds between G^ and S^. 



Let us now apply this general formalism to the calculation of the quantum corrections 
discussed in Sect. 3.2.1. We then need the following contour Green's functions: 

(a) The two-point function G'^^j^" {x ^ y) = {5Al{y)SA'^{x)) of the soft fields induced by 
their coupling to the semi-fast gluons. This is given by eq. ([4.31|) with replaced by 

SE<,ix,y) = {S.My)S.Ux)) , (4.32) 

where SJ^ is the quantum color current in eq. ( |3.14| ). To lowest order in , the soft field 
propagators Gji and Ga in eq. ( [4. 311) can be computed in the mean field approximation; 
that is, one can neglect the quantum self-energy STi'^ when solving eq. ( |4.3CI| ). With these 
approximations, eq. ([4.31|) reduces to eq. ( p.l5| ), which has been used in Sect. 3.2.1. 

(6) The background field propagator iG(^{x,y) = {Tc a^{x)a^ {y)) of the semi-fast 
gluons in the Gaussian, or mean field, approximation (cf. eq. ( p.82|) ). This is given by 
eq. ( |4.28| ) where the quantum self-energy 5S^j, is neglected; that is, 

V,,{x)G''c\x,y) + V / d^zlf{x,z)G-c\z,y) = 5l5c{x,y) . (4.33) 

Since the tree-level self-energy ( 4.25 ) is non-local, the integration over z^ in this equation 
runs along the whole contour G. Thus, even for real time arguments x"*" and y"*", eq. ( |4.33| ) 
is still sensitive to the complex side C_ of the contour. This is to be contrasted to 
a local field theory, where the contour structure becomes irrelevant in the mean field 
approximation ||4^ . 

However, a closer inspection reveals that, even for eq. i WM ). the contour structure is 
not essential, in the sense that one can ignore C_ when computing the Feynman propa- 
gator (i.e., the propagator G'^{x,y) with both time arguments on C+). Indeed, for x~^, 
y^ E G^, the contribution of C_ to eq. ( |4.33| ) reads, schematically, 



d'^zU^{x,z)Gc{z,y) = - J d^zU<{x,z)G>{z,y) oc 

J dz^ G>{z^ -y^,x,y) = G>{p~ = 0,x,y) = 0, (4.34) 



ic. 

oc 



where we have used eqs. ( [4.25| )-( ^.26| ), the time- homogeneity of the problem (i.e., the fact 
that II"^ is independent of time, while G^ depends only upon the relative time z^ — y^), 
and the fact that G^{p~ = 0) = (since, by definition, this has support only for in 
the strip (^M) )- 

To conclude, the background field Feynman propagator iG'^'^{x,y) = (T a^{x)a'^{y)) 
can be obtained by solving the following equation 



{x)G-'\x,y) + S,^Jd'zU{x,z)G-\z,y) = 6i6^^\x - y) , (4.35) 
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where all the time variables are real, and n(x, z) is the restriction of Il'~''{x, z), eq. ( |4.25| ), 
to time arguments on C_|_. Then, the corresponding Wightman functions (which enter a 
and X ; see, e.g., eq. ( 4.17|) ) can be extracted from the Feynman propagator by considering 
appropriate time orderings, as we did for the free propagator in eqs. (|4.19|) - (|4.21|) . 

Note that any explicit use of the contour has been avoided in this way. This has 
been possible because, first, of the static character of the background, and, second, of 
the simple nature of the present, mean field, approximations. For more general problems, 
where higher-order quantum effects should be computed in the presence of inhomogeneities 
in time, a complete use of the contour would be generally unavoidable. 



4.3 Feynman rules for a and x 



As explained in the previous subsection, the use of the contour techniques can be avoided 
for the calculation of a and x- It is then sufficient to consider the restriction of the 
quantum theory to the real time axis, as defined by the following action (compare to 
eq. (U)): S = Sym + Sw, with 



w 



d'STi [p{x)Woo,-oo[A-]{x)] 



and the real-time Wilson line: 



\A-]{x) = T exp 



(4.36) 



(4.37) 



This is the original version of the quantum MV model proposed in Ref. |Tl[]. In particular, 
the tree-level self-energy generated by this action reads: 



W\x,y) 



w 



5A-{x)5A, (y) 
gp'^'ix) 5^'\x - y) {x^\FV ^ly-') 



r^>,(f)e(x+-y+)(5(3)(f-y) 



(4.38) 



where p"* = — i/^'^^p'^, and e{x) = 6{x) — 9{—x). Eq. ( [4.38| ) coincides, as it should, with 
the restriction of the contour self-energy 11'^, eq. ([4.25|) , to the real time axis. Note the 
emergence of the principal value prescription in l/p~ (this is defined as in eq. ( p.85| )). 

Since we are mainly interested in the strong field regime, where ^ 1/g and p ~ 1/g, 
it is convenient for power counting to perform a rescaling of the background fields and 
the color sources: 



At: 



9 



Pa 



9 



Pa, 



Spa 



9 



Spa 



(4.39) 



In the saturation regime, the new tree-level quantities A'^ and p are of order 1, while the 
new quantum corrections x ^"^^ — which are defined in terms of the rescaled charge 
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fluctuations 5p by the same equations as before (i.e., ( |3.17] ) and ( |3.23| )) — are of order 
Ug. After this rescaling, the coupling constant g completely drops out from the classical 
analysis in Sect. 2. On the other hand, we prefer not to rescale the quantum fluctuations 
; thus, their propagator G^^{x^y)[A, p] remains of order one, and all the factors of g 
can be explicitly read off the vertices (see below). 



The quantum charge fluctuations 5pa{x) are defined by eq. (|3.14|) with p = — . After 
the rescaling ([4.39|) , this yields: 

5p„(x) = 5p«(x) + 5pf (x), (4.40) 

where 5p^^^ is linear in a^, and dp^"^^ is quadratic: 

5p«(x) = -2z(?^+;(f)a-(a;) + 

+gp^'{x) J dy+{x^\PV-^\y+)a'^-{y+,x), (4.41) 

5p^^\x) = g'r'"^[d-'a1ixMix) 

2 

- j;^p\x) I dy+a-\y+,x) j dz+a-\z^,x) 

X |^(z+ - y+)d{y+ - x+)Tr (t^t^T'^T^) 
+ e{x^ - z^)e{z+ - y+)Tr {T'^T^T^T'^) 
+ e(z+ - x+)e{x+ - y+)Tr {t^T'^t^T^) | , (4.42) 

where all the factors of g are now explicit. In the right hand sides of these equations, the 
terms involving a* come from the three-gluon vertex in Sym-, while the terms involving 
a~ come from the two- and three-point vertices in Sw- 

In eqs. ([4.41|) - ([4.42|) it is understood that only the soft modes (with k'^ ^ ftA"*") are 
kept in the products of fields. This has been used to simplify the first contribution to 
(5p(^) by writing: 

(9+^L) a'" - d^a'' ^ 2 J-,+>- . (4.43) 

Thus, the above expressions can be used only for calculations to LLA. 

The quantum corrections x and a are given by eq. ( p.l7| ) (with x"*" = y^) and ( |3.23| ), 
where the average over the semi-fast fluctuations is defined as in eqs. (|3.82|) - (|3.83|) . 
This generates Wightman functions of the fields (see, e.g., eq. ( [4.17| )), which are related 
to the corresponding time-ordered, or Feynman, propagator as explained in Sect. 4.2. In 
particular, the equal time limit of a two-point function is obtained as: 

{alix^Xix"-)) = G<r(x+ = 1/+) = ^GZ{y'' = x+ + e). (4.44) 

In momentum space, this instructs us to perform the integration over p~ by closing the 
contour in the upper half of the complex p~ plane: 

G'<(a;+,x+) = I j ^e'P"^G(p-). (4.45) 
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a) 



b) 




c) 



d) 



Figure 5: Feynman diagrams for o"i (a) and 02 (b,c,d). The wavy line with a blob denotes 
the background field propagator of the semi-fast gluons; the continuous line represents 
the source p; the precise vertices can be read off eq. (|4.42|) . 



In practice, this prescription will play a role only for the vacuum piece Go (i-e., the limit 
of G as p — >■ 0, which to the order of interest coincides with the free propagator Go; cf. 
Sect. 3.4). Indeed, within the present approximations, the induced piece bG = G — Gq is 
continuous at = ■, so its equal-time limit is unambiguous. 

We are now prepared to express a and x terms of the Feynman propagator G^^(x, y) 
of the semi-fast gluons. For a = (5p) = (^p'-^^), one obtains (with the "approximately 
equal" sign denoting an equality which holds to LLA) : 



asln^cra(x_L) ^ ^ dx Tr(T"o-(f)), 



a{x) 



g'd^G^' + ig'p{x)ix 



x=y 

= (Ti(x) + a2(x). 



id' 



G- 



id~ 



X 



(4.46) 



In writing o"2 as above, we have used compact but formal notations for the second, non- 
local, contribution to 5p^^\ eq. ( [4.42| ). A pictorial representation of the two pieces of a is 
given in Fig. ^. 

Concerning x, we note that to order this involves 5p^^\ but not 5p^'^\ [This is clear 
by counting powers of g according to eqs. ( [4.41| ) and ( [4.42] ).] Therefore (with = x~^ + e): 



dx / dy Xab{x,y), 



Xab{x,y) 



9 



liT^^d' + p 



id 



id' 



(4.47) 



where we have used also the symmetry property 

1 , , 



(x|PV 



id~ 



id' 



\x). 



(4.48) 
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c) d) 

Figure 6: Feynman diagrams for the four contributions to x given in eq. ( [4.49| ). 



This further yields, in matrix notations (the PV prescription in 1/p is imphcit) : 



1 



+1 

X 



px {A 



id' 



G' 



+ ipx{x\ 



id' 



-G- 



)Py 
1 

id~ 



Py 



(4.49) 



Diagramatically, all the above contributions to x ^i-re represented by tree-like Feynman 
graphs (no loops), with vertices proportional to p or (see Fig. Thus, like the 



tree-level source p itself, x{^i v) is localized near the LC, at ^ x , ?/ 



< 



1/A^ 



It is now easy to verify, by using eqs. ( [4.40| ) and ( [4.42| ) and straightforward power 
counting, that all the n-point correlators of 5p beyond n = 2 are of higher order in as- 
Indeed, each factor of in these equations is accompanied by a power of g. Then, for 
instance, the three-point function {6p6p6p) — which involves, at least, four fields — 
is of order a^. This explains why, to the present accuracy, we can restrict ourselves to the 
one- and two-point functions in eqs. ([4.46|) and (|4.49| ). The remaining part of this paper, 
and also Paper II, will be devoted to their explicit evaluation. 
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5 The weak source approximation and BFKL 



As a simple illustration of the formalism developed in the previous sections, let us explic- 
itly compute here some of the contributions to a and x the weak source (p <C 1), or 



linear response, approximation. We already know from Ref. that, in this approxima- 
tion, eq. (|3.56|) reduces to the BFKL equation. Our intention here is not to repeat the 
derivation in Ref. [|T^, but rather to illustrate in this simple setting some of the general 
arguments which have appeared in the previous discussion. These include the longitudinal 
structure of the induced source, the equivalence between the strip restrictions on p"*" and 
on and the sensitivity to the ie prescription in the LC-gauge propagator. 

To this aim, we shall compute just one contribution to a and one contribution to Xi 
namely (cf. eqs. ( |4.46| ) and ( [4.49| )): 



(5.1) 



These quantities will be evaluated here to lowest order in p (that is, to linear order in 
the case of ex, and to quadratic order for x), and to LLA. In this approximation, there is 
no distinction between the LC-gauge source p and the COV-gauge source p, so we shall 
use the notation p to refer to any of them. Also, to simplify the notations, we shall not 
distinguish below between a and a, neither between x ^"^^ X- 



5.1 The weak source hmit of a 



Consider cti first: because of the color trace in eq. ( [S.!] ), the vacuum piece Gq of the 
propagator drops out from this equation, which can thus be rewritten as 



-g'd^Tv{T^5G'\x,y))l^^ = |-(a+ - d^)Tr[T'^5G'\x,y 



x=y 



(5.2) 



where the second equality follows from the symmetry property SG^^i^{x,y) = 6Gll^{y,x). 

A contribution linear in p to 56*** can be generated either by a direct insertion of p, 
which couples to the fields a~ via the vertices in Sw (cf. eqs. (^4.36|) - (^4.37|) ), or by an 
insertion of the mean field A^, which couples to a' via the usual vertices in Sym- However, 
the insertion of does not contribute to cxi to LLA: indeed, this produces the standard 
one-loop polarization tensor of QCD, which has no logarithmic enhancement. On the 
other hand, such a logarithm is generated by the insertion of p, as we show now. 



The corresponding contribution to (Xi is displayed in Fig. ||.a, and reads 



ai[x) 



-g'l I d^yGi-ix-z)Iliz,u)d^G^'iu-y) 



x=y 



(5.3) 



where n(z, u) is given by eq. ([4.38|) , and Gq and Gq * are the free propagators in eq. (|3.84|) . 
Below, we shall compute this quantity in two different ways: 
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a) First, we shall directly compute the two-dimensional charge density <7a{x±) = J dx~(Ta{x) 
for which the logaritmic enhancement is manifest (cf. eq. (|4.46| )). We shall verify that 
the result is independent upon the choice of the strip restriction, and also upon the ie 
prescription in the LC-gauge propagator. 

b) By using a strip restriction on p~ , we shall recalculate ai{x) in the x~ representation, 
and thus exhibit its longitudinal structure. 

a) We first rewrite eq. (|5^ ) in momentum space, as follows (with q = p — k): 



PV — p(g-) 
P 



k+G-,\p-,k). 



Here, k^ could be as well replaced by or (A;+ +p+)/2 : this can be seen by using the 
symmetry of the integrand under the simultaneous exchange k ^ —p and p" —p~, 
or, alternatively, by using directly the second, more symmetrical, representation for cxi in 
eq. (|5.2| ). By also using eq. ( p.84|) , we can write: 

- ^(;^ + I^)go(p)GoW. (5.4) 

The integration over x~ then sets k~^ = P~^, so we are left with the following, manifestly 
symmetrical, integrations over p~^ and p~ : 

^PV^ /^PV- ^ ^ (5.5) 

27r p~ J 2n p+ 2p+p~ — p\ + ie 2p+p~ — k\ + ie 

Note the effective PV- prescription in this is a consequence of the symmetry proper- 

ties of the integrand in eq. (^.4| ), and is largerly insensitive to the original ie prescription 
in the LC-gauge propagator ( p.84|) . This same effective prescription would have been 



generated by starting, e.g., with the advanced prescription of Refs. [^, 0, or with the 
principal value prescription. Thus, the present calculation of ai is independent of the 
axial prescription. 

The integral in eq. (|5.5| ) must be supplemented with a strip restriction either on p"*" 
(cf. eq. (|3.8|)), or on p~ (cf. eq. ( |3.8U| )), with the two restrictions related by eq. ( |3.11| ). 



Given the symmetry of the integrand under p"*" ^ p~, and the fact that the integral is 
only logarithmically sensitive to its cutoff, it is clear that both restrictions will lead to the 
same result. For instance, with a strip restriction on p^, the unrestricted integral over p~ 
can be performed by contour techniques, to yield: 

PV ^ ^ ^ - '<P^^ . (5.6) 

27r p^ 2p+p^ — p\ + ie 2p+p^ — k\ + ie 2p\k\ 

Then, the restricted integration over p"*" generates the expected logarithmic enhancement: 



strip 27r p+ W-A+ JhK+ / 27r p+ vr 
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One thus obtains: 



g'^Nc 1 r (fp± f (fk 
^hJ (2^ J (27) 



Pi. ■ fcj 



in agreement with Ref. [JIOl- In this equation, 

Pa(g±) =Pa(g^ = 0,gx) = dx' d^x±e''^^-''^ pa{x',x±) . 



{5.t 



(5.9) 



After simple manipulations, eq. ( ^.8|) can be finally rewritten as: 



<(gx) = -2Area,ln^p,(gx) / 



2plip±-q±)^ pi 



(5.10) 



b) The longitudinal structure of eq. (^.3|) is contained in (cf. eq. ( [4.38|) ) 

' dz- Gl-{x- - z-)p{z-)dtG,\z- ~x-) ^ pGl-{x-)dtG,\-x-), (5.11) 



where we have used p{z ,2;^) ~ 5{z )p{z±), which is appropriate since p{z ) is localized 
at ^ 2;^ ^ l/A"*", while <7{x~) is rather needed at large x' ^ 1/6A+. 

With the strip restriction ( p.80| ) on p~ , the integrations over p'^ and are unre- 
stricted, so we can use the standard expressions for the propagators Gq~ and Gq * in the 
x~ representation]^ : 



dp'^ 



\ 9{x-) - \e{x-)e{p-) 
PI I ^ 



p+ + ie 2p+p ~ p\ + ie 

i{-x-)e{-p- 



(5.12) 



— I — —X 

e 2p- 



d+G^'i-x ,p ,k^) 
This yields then: 



2p- 



[e{-x~)e{p-) -e{x-)e{-p-)}e^"' . (5.13) 



Gl{x-)d:G,\-x-) = -.^(x-)\^^^e4 

zp p^ 



(5.14) 



which implies that ai{x ) has support at 1/A+ < x ^ 1/6A+. Indeed, the typical range 
of x~ is fixed by the exponential e^^^ff^ , with p^^j = k\/2p~] since p~ is restricted to the 
strip ( |3.80| ), p^jj is constrained by 6A+ < Ip^IjjI ^ A+, hencefrom the support of ai( 
alluded to before. 



(X 



^^By contrast, it would be more difficult to study the a; structure by using a strip restriction on 
cf. the discussion in Sect. 3.4. 
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In fact, in this case it is easy to pursue the calculation and obtain ai{x^) in explicit 
form. After inserting eq. ( |5.14|) into eq. (|5.3| ), one has to compute the restricted integral 
over p~ . This yields: 

'strip 2tx 2(p-)2 ^ ~ kl 27iix ' ^ ■ 

where we have replaced k']_/2Ar A"*" to LLA (cf. eq. (|3.11|) ). The final expression for 
cri(x~) reads then (compare to eq. (|5.8|)) : 



with the "form factor" 

Fix-) = e(x-)- , (5.17) 

X 

which is indeed localized at 1/A+ ^ x^ 1/6A+, and yields the expected logarithmic 
enhancement after the integration over x~: 



I 



dx-F{x-) = In^. (5.18) 



Of course, the specific longitudinal structure in eqs. (|5.16| )-( p.l7| ) is directly related to 
our choice of the retarded prescription in the LC-gauge propagator ( |3.84|) . If the advanced 
prescription were chosen instead, the corresponding ai{x~) would be localized at negative 
x~ , while with a PV-prescription it would have support at both a;~ > and a;~ < 
(with 1/A+ ^ \x^\ ^ l/bA~^ in all cases). Nevertheless, the integrated charge density 
o'i{x±) = / dx^cri{x~ ,x±) — which, we recall, is the relevant quantity to LLA — , would 
come up the same with any prescription, as already demonstrated by calculation (a). 

To linear oder in p, the second contribution to the induced charge ( [4.46| ) reads 

(72(x) = -ig'^p{x) {x\^Gq'^\x), (5.19) 

and involves the prescription-independent component Gq~ of the free LC-gauge propaga- 
tor. To evaluate eq. ( |5.19|) , one has to recall the precise temporal and color structure of 
eq. (|4.42|) . After some simple algebra, one obtains |T0( 

a-^{q^) = -2N,aM^pa{qA_) j (5-20) 

which precisely cancels the tadpole piece 1/px ^^e integrand of eq. (|5.10|) . Thus, finally, 

^ < + = -iV.„.lnip.(,,)/i^-jj-^, (5.21) 

where the upper script on a'-^^ is to recall that this is the result to linear oder in p. This 
result holds independently of the axial prescription. 
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5.2 The weak source limit of x 



To lowest, quadratic, order in p, we have (see Fig. |l[a) : 

Xt{x,y) = WJ^:::{x){x\G'^\y)J^:,\y), (5.22) 

where = + e, and is related to p by the linear-response, or Abelian, approxi- 
mation (cf. Sect. 2.2): ^ -{d'/WDp. Eq. ( |0^ involves: 

^ ' i 27r i 27r i (27r)2 2p+j9--pi+ie ^ ^ 

which we shall first evaluate with a strip restriction on p"*". Then, e.g., ^ 1 (since 

typically \p'^\ <C A"*", and a;~ ^ 1/A"^), so we can neglect x~ and ?/~ in eq. ( |5.23| ). By 
using the convergency factor e*^ the integral over p~ can be computed by closing the 
contour in the upper half plane. This yields: 



Xi{x,y) = -Ag^J^^\x 



dp'^ 9{—p^) f d^p±_ 



strip 27r 2p+ J {2ti 



e 



ipx-(xx-yx) 



= ^\xi{l/h)5^''\x^-yA_)T^\x)T^\y). (5.24) 
vr 

It is a little bit trickier, but at the same time instructive, to show that the same result 
can be obtained by using a strip restriction on p~ . Then, the convergency factor plays no 
role, but the (restricted) integration over p~ can still be performed in the first place, by 
using: 

'-f « ^ e.«,(p-). (5.25) 

strip 27r 2p+p —p^j_ + ie Ap+ 

which holds to LLA. Here, Qstrip{p^) is a step function which implements the strip restric- 
tion (|3.8|): Qstrip{p^) = 1 for feA"*" <C -C A+, and Qstrip{p'^) = otherwise. Thus, to 
LLA, the restricted integration over p~ generates an effective strip restriction on p~^, with 
the two restrictions related by eq. ( |3.11D . This is so since, to this accuracy, the integral 
in eq. ( |5.25| ) is saturated by on-shell modes. It is now easy to verify that the ensuing 
integral over p+ leads to eq. (|5.24| ) again. 

The other components of x in sq. ( [4.49 ) can be similarly computed, with the following 
final result (in matrix notations where, e.g., p^ = Pab{x)) [p!0[] : 

X'^\x^y) = 4«.ln(l/6)/|^^^^ 

X {p.py + ip.{T^Y)y - i{T^y).Py + p\Tt'T^') . (5.26) 

The result (|5.26|) is independent of the axial prescription in the LC-gauge propagator, as 
can be easily verified by repeating the calculations with various prescriptions. 
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5.3 The BFKL equation 



In the present weak field approximation, the only non-trivial evolution equation is that 
for the 2-point function {pp)t, which reads (cf. eq. ( p.56|) , which we also integrate over 
x~ and y~): 

^{pa{x±)pbiy±))^ = as{cr^a\x^)pb{y±) + Paix±)ai^\y±_) + Xabi^±,y±))^- (5-27) 

It is interesting to consider the diagonal element of this equation in momentum space. 
Indeed, this is the equation satisfied by 

Pa{k±)pa{-k^))^ ^ kl{\K-'{k^)\')^ = V{x,kl), (5.28) 

which, according to eq. ( p.39|) , is proportional to the density of gluons with transverse 
momentum k± and rapidity r = ln(l/x) (i.e., longitudinal momentum k~^ = xP^). This 
quantity is often referred to as the unintegrated gluon distribution. Note that the simple 
relation (|5.28| ) between ip{x, k"]^ and (pp)r holds only in the weak source limit, where the 



electric field and the color source are linearly related: —i[k^/k\)p. 

By combining eqs. ( |5.27|) - (|5.28| ) with the previous expressions for eq. { ^ ), and 



X^'^\ eq- ( |5.26| ), one finally obtains the following evolution equation for ip{x, k'jj 
dip{x,kl) asNc f ,2 kl 



X- 



f d^pi_ ^ -^{x,p\) 

J pi{k±-p±y 

9 2 fd'p^^Tr^ T^vi^,kl), (5.29) 

27r^ J Pj_{k± — p±) 



dx TT^ J P±ik± — p± 

a,Nr, r ,n k^ 



which coincides, as anticipated, with the BFKL equation [^. The first term in the r.h.s., 
which here is generated by x^'^^ ? is the real BFKL kernel, while the second term, coming 
from a^^^ , is the corresponding virtual kernel. 



Thus, outside the saturation regime, the general evolution equation (|3.53|) reduces to 



the BFKL equation, as necessary on physical grounds. This is a very non-trivial check of 
the quantum structure of the effective theory under consideration (i.e., eq. (|3.1j ) with the 
action ( ^4.8| )). Indeed, unlike the classical MV model, which is sensitive only to the eikonal 
vertex / d'^xpA~ (the one-point vertex in the action Sw), the quantum corrections a and 
X are generated by the 2-point and 3-point vertices in Sw- These are also the vertices 
which enter the calculation of a and x the general non-linear case, to be described in 
Paper II. 

As we have explicitly verified here, the BFKL equation is obtained independently of 
the gauge-fixing prescription in the LC-gauge propagator. This is as expected since, in this 
approximation, the 2-point function {pp)r of the color charge coincides with the physical 
gluon density, cf. eq. ( ^.28|) . In the presence of non-linear effects, this relation may not 



hold anylonger. Then, the charge correlators are not physical observables by themselves. 
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and need not be gauge- fixing independent. And, indeed, we shall find in Paper II that 
the non-linear effects in a and x ^'^^ sensitive to the axial prescription. 

To anticipate, the simplest results will be obtained by using the retarded, or the ad- 
vanced, prescriptions discussed in relation with eq. (|3.84|) . With the retarded prescription, 
both the classical mean field eq. ( p.31| ), and the induced source cr, eq. ( |5.16| ), sit at 



> 0, or z < t ; that is, the soft gluons fields are behind their source, the hadron 
(which is located aX z = t). It is intuitively plausible that within this picture there are no 
initial state gluon interactions. This is the counterpart of the conclusion of Mueller and 
Kovchegov |l5l who used an advanced prescription — ie), and showed that all the 



final state gluon interactions disappeared. These two prescriptions are simple in that one 
can either put gluon interactions in the final or initial state. Other gauge conditions such 
as Leibbrandt-Mandelstam or principal value do not have this simple feature, and lead to 



other complications for computations as well |32 



6 The background field gluon propagator 

In this section, we construct the background field propagator iG^^{x,y) = {T a^{x)a''{y)) 
of the semi-fast gluons, which is the key ingredient for the the calculation of a and x 



to be presented in Paper II. (See also Refs. [^, ^ ^ for previous related work.) The 



propagator will be obtained by solving eq. ( |4.35| ) in the LC gauge = 0, and in the sub- 



space of fields with energies p~ restricted to the strip (|3.8CI|) . The latter restriction causes 
no difficulty, since the problem being homogeneous in time (recall that the background 
field is static), we can construct the propagator for any fixed p^. Concerning the gauge 
condition, we still need to fix the residual gauge freedom in the LC gauge; this will be 
done by giving a prescription for the "axial" pole at p"*" = 0, in agreement with eq. ( |3.84| ). 

For more clarity, we shall first briefiy review the corresponding construction for the 
scalar propagator, whose definitive form has been given by Hebecker and Weigert ||37|| . 



Then, in Sects. 6.2 and 6.3, we shall construct the gluon propagator in two steps: First, 
from the scalar propagator, we shall derive the gluon propagator in the temporal gauge 
a" = 0. Then, the propagator in the LC gauge will be obtained via a gauge rotation. 
The reasons for this particular strategy should become clear in a moment. 



6.1 The background field scalar propagator 

The background field scalar propagator G{x,y)[A] is defined as the solution to 

-VlG{x,y) = 6^'\x-y), (6.1) 
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whereQP^ = V+V' +V-V+ - Vl = 29+^- -Pi , ^"^ = d^'-iA^', and = S^^'A' is the 
LC-gauge color field generated by p and constructed in Sect. 2.3. To solve this equation, 
it is convenient to perform a rotation to the covariant gauge where the background field 
has just a plus component: A'^ ^ = 5^^a (cf. eq. (|2.21|) ). We thus write: 

G{x,y) = U{x)G{x,y)U\y), (6.2) 

where G{x,y) obeys the same eq. ( |6.1| ), but with = 2d^V^ — V^ and "D^ = — ia{x). 

In order to solve this equation non-perturbatively in a, it is important to recall the 
separation of scales in the problem: the semi-fast quantum fluctuations whose propagator 
we are about to compute have relatively small longitudinal momenta (typically, in the strip 
( p.8| )), and therefore a resolution scale Ax~ which is much larger than the longitudinal 
extent ~ l/A"*" of the source p. When "seen" by these fluctuations, the source p and the 
associated field a are effectively delta functions at x~ = 0. We are thus led to solve the 
following, singular, equation: 

{-2d-[d+ -i6{x-)a{x^)) +Vl}G{x,y) = 6^^\x - y) , (6.3) 

whose precise meaning is as follows: for x~ either positive or negative, G{x, y) obeys the 
free equation 

-dlG{x,y) = 5^^\x-y) for a;" ^ 0, (6.4) 

but it is discontinuous at = (and similarly at = 0). To evaluate its discontinuity, 
it is easier to argue in the LC gauge: the transverse background field A^ is effectively a 
step function (cf. eq. ( p.31| )), so it is discontinuous but finite at x~ = 0. Then, eq. (|0| ) 
implies that G{x, y) must be continuous at a;~ = : 

Gix,y) U-^o+= G{x,y) \x-^o-, (6.5) 



and similarly at y =0. After the gauge rotation ( |6.2| ), this provides the sought for 
boundary conditions for G{x, y) : 



=0") 



G{x,y) \^~=o+ = V\x±)G{x,y) 
G{x,y)\y-=o+ = G{x,y) \y-=o- V{y±) , (6.6) 

with V^{x^) defined in eqs. (^[^-(^[SUp. 

To summarize, to construct the scalar Green's function, one must solve the free equa- 
tion ( |6.4|) subject to the boundary conditions ( |6^ . The solution reads : 

G{x,y) = Go{x,y) -2 J d'^zGo{x,z)S{z~)d~Go{z,y) 

X {eix-)ei-y-)iV\z^) - 1) - ei-x-)eiy-)iViz^) - l)} , (6.7) 

^^Recall th at th e gauge couphng constant has been absorbed into the normahzation of the classical 
field, cf. eq. (^!39|) . 
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Figure 7: The propagation of the scalar particle in the delta potential at = 0: a) a 
crossing trajectory {y~ < and x~ > 0; b) a non-crossing trajectory (both x~ and y~ are 
positive); c) a reflexion on the potential; this last process is actually forbidden by energy 
conservation. 

where Go is the free scalar propagator: 

To verify that this is the correct solution indeed, note that: a) In the absence of the 
background field, V = V'^ = 1, and G reduces to Gq, as it should, b) When acting on 
eq. ( |6.7| ) with d^, one gets zero everywhere except at = when operating on the 
left, or = when operating on the right. Therefore G solve the equations of motion 
everywhere except where the source sits, c) By using 

Goix) = -^ J {9ix-)9ip-)-9i-x-)9i-p-)} J ^e^^x --^.-^ (g g) 

one deduces that 

\im^ Goix-z) U-=o=-tS^'\x^-y^) J ^ ' ^ {^(x-)^(p-) - 0(-a;-)0(-p-)} 

(6.10) 

Inserting this into eq. ( |6.7|) shows that G satisfies the correct boundary conditions across 
the source. 

After a Fourier transform to the energy representation, eq. ( |6.7| ) can be rewritten as: 

G{x,y,p-) = Go{x-y,p~){9{x-)9{y-) + 9{-x-)9{-y-)} 
+ 2ip' J d^zGo{x-z,p')5{z')Go{z-y,p') 

X {9{x-)9{-y~)V^{z^) - 9{-x-)9{y-)V{z^)} . (6.11) 
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This result is quite intuitive. There are two types of terms (see Figs, and b): A 
crossing term, where the endpoints x~ and y~ are on opposite sides with respect to the 
plane x~ = (where the potential a is located), and a noncrossing term, where both 
x~ and y~ are on the same side of the potential. In the latter case, the scalar particles 
propagates freely from y to x . There is no reflection on the potential, i.e., no process like 
the one depicted in Fig. ^.c ; this is forbidden by energy conservation, since the potential 
is static. Indeed, for on-shell excitations, and p~ have the same sign; thus, a particle 
with positive energy propagates forward in x~ , as obvious on eq. ( |6.9| ), and this cannot 
be changed by its scattering off the static potential, which preserves p~ (see Fig. 0.b). In 
particular, this shows that one can effectively replace 

9{x-)9{-y-) ^ 9{p-), 9{-x-)9{y~) 9{-p-), (6.12) 

without changing the crossing term in eq. ( |6.11| ). 

After also performing the rotation ( |6.2|) to the LC gauge for the background field, we 
finally derive the following expression for the background field scalar propagator: 

G{x,y,p~) = Go{x-y,p-){9{x-)9{y-)V{x^)V^{y^)+9{-x~)9{-y-)} 
+ 2^p-Jd^,Goix-^,p-)6i^~)G.i^-y,p^) 

X {9{p-)V{x^)V^{z^) - 9{-p-)V{z^)V^{y^)}. (6.13) 

It can be easily verified that this function is continuous at both x^ = and y~ = 0. 



6.2 The gluon propagator in the temporal gauge 

To construct the propagator in the light cone gauge a"*" = 0, it is convenient to perform an 
excursion through the temporal gauge d~ = 0, that is, to compute first the propagator in 
a mixed gauge where the background field is in the LC-gauge0 = 0, but the quantum 
field is in the gauge d" = 0, and then rotate the quantum field to the LC gauge a+ = 0. 
This last rotation will be constructed in Sect. 6.3 below. Here, we shall rather focus on 
the construction of the propagator in the temporal gauge. 

Let us first explain the usefulness of such a detour. The reason is technical: by 
inspection of the coefficients in eq. ( |4.35|) for G'^'^, that is (cf. eqs. ( [4.24| ) and ( |4.38| )), 



G;l{x,y)^{g,,V'-V^V,-2zJ^^,)/^\x-y) - 6,^6,^- p{x) e{x+ - y+) 6^'\x - y), 

(6.14) 

one finds that some of the coefficients which count for the LC gauge = 0, namely 
Gz^, GZi and G^}, become singular in the limit in which we treat p as a delta function, 

^°Unless otherwise specified, the background field will always be in the LC-gauge in what follows. 
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p{x) 6{x^)p{x±). (These coefficients involve either p, or JF+*, which is itself singular 
in this limit, cf. eq. (|2.32|) .) On the other hand, there is no such a singularity in any 
of the coefficients relevant for the gauge d~ = 0. Thus, in the latter gauge, the vector 
propagator G'^'^ (with p,^ = +,1 or 2) is continuous at a;~ = and y~ = 0, and can 
be immediately obtained from the scalar propagator discussed in the previous subsection. 
This construction is described now (see also Refs. [^, |3^ for more details): 

To simplify notations, in the remaining part of this subsection we shall denote the 
quantum fields in the temporal gauge simply as (thus, a~ = 0). The easiest way to 
obtain G'^'^ is to use its definition in terms of a path integral (the analog of eq. (|3.82|) ): 

iGZ{x,y) = Z~'Jva6{a')a^,{x)al{y)e''^^^''^\ 

So[A,a] = J d^xi^^a\-V^)a' + ^{d-a^ -V'a'f^ (6-15) 

where the second line follows by evaluating the Gaussian action ( p.83| ) in the temporal 
gauge a~ = 0. Note that this new action does not depend upon p explicitly, but only 
implicitly, via the background field A^. This action can be brought to a quadratic form 
via the redefiniton 

a+ = a+ + -^W. (6.16) 

There is no ambiguity associated with the pole at p~ = in this transformation since, by 
definition, the LC energy is restricted to the strip (|3.80|) . We thus obtain: 

S[A,d+,a'] = J d^xl^^a'{-V^)a' + ^{d'a+^Y (^-l^) 

which immediately implies that the propagator G*-' of the transverse fields is the same as 
the scalar propagator: 

G'\x,y) = 6'^Gix,y) (6.18) 
with G given by eq. ( |6.13| ). Also, 

(T d^xra^y)) = / ^ e--(^-) ^ , (6.19) 

while, e.g., (T a^(a;)a*(y)) = 0. Thus, the mixed propagators read 

G+\x, y) = -^ViG{x, y), G^+(x, y) = G(x, y) ^ ^ (6.20) 

where T)^ = d + igA, with the derivative d acting on the functions on its left. The ++ 
propagator includes also the contact term (|6.19| ), and reads: 

G++(x,i/) = ^VlGix,y)Vl^^ + l-l-5i^){x-y). (6.21) 

^x ^y 
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These propagators are conveniently rewritten in the energy representation: 

G^\x,y,p-) = ±VlG{x,y,p-), &^{x,y,p-) = - ± G{x,y,p-)Vl^ , 
p p " 

G++ix,y,p-) = {viG{x,y,p-)Vl^ + 5^^\x - y)] . (6.22) 

We stress again that, in our p~ renormahzation group scheme, we remove exphcitly the 
support of the quantum fields at p^ = 0, so that the operator l/p~ is not singular. 



6.3 The gluon propagator in the LC gauge 

The vector propagator G^'^ is defined by eq. ( p.82|) where both the background field and 
the quantum field are in the light cone gauge: A'^ = a"*" = 0. It is however straightforward 
to perform a change of variables in the form of a gauge rotation, so as to reexpress G'^'^ 
as an integral over quantum fields in the temporal gauge d~ = 0, and thus relate it to the 
propagator G^'^ in the previous subsection. 

Specifically, the gauge rotation d^ reads: 

A^^af" = A{x)(A^ + a'' + id'')A\x), (6.23) 



with (the initial time Xq is arbitrary): 

A^(a;) = Texpji^ dz+ a' {z+ , x)'^ . (6.24) 

Since A~ = 0, this transformation insures that d~ = 0, as it should. For consistency with 
eqs. ( |3.82| )-( p.83D , where the small fiuctuations are treated in the Gaussian approximation, 
the gauge rotation ( |6.24|) must be expanded to linear order in a~ : 

A^(x) ^ 1 + iA(x), d'A{x) = a~{x). (6.25) 

This implies, in particular, that the field A{x) carries the same energy p~ as the quantum 
field a~ , i.e., a p~ restricted to the strip ( p.SOp . As a consequence, 

A(x+) = / — e'P""^A(p-) (6.26) 

J strip 

vanishes asymptotically, A(a;"'") — as |a;^| oo, so the infinitesimal gauge transfor- 
mation ( |6.25| ) satisfies the periodicity condition /\{x^ —>■ oo) = A{x~^ —>■ — oo) = 1, as 
required for this to be a symmetry of the action ( [4.36| ). 

To linear order in the small fiuctuations, eq. ( |6.23| ) reduces to 

d^ = - [V^", A] . (6.27) 
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This is a symmetry transformation of the full action ( ^.361 ), and also of the quadratic 
action (|3.83| ). [To verify the latter property explicitly, use the gauge invariance of the full 
action, together with the fact that A'^ is a solution to the classical equations of motion 
6S/6A = 0, to deduce that {6'^S/6A^'6A'')aV''A = for arbitrary A .] 

Eq. ( |6.27| ) with A defined by eq. ( |6.25| ) gives us the change of variables to be performed 
in eq. ( p.82|) . To implement this, we first rewrite eq. ( p.82|) as 



JVaJVA6{d-A 



JVaJVA5{d-A 



iSo [A,a] 



(6.28) 



where the functional integral over A is just to relabel a a.s d A. (This transformation is 
non-singular since p~ never vanishes.) Then, we change variables according to eq. ( |6.27| ): 



+ [P^ A] 



(6.29) 



and perform the integral over A with the help of the delta function S{a~^) = 6{d~^ + d'^A). 
This finally yields: 



J Vd5{d-) {d^' + V^'A)^{d'' + T>''A)^e'^°^-^'^^ 



JVd6{d~) e^'5o[Aa] 



(6.30) 



where it is now understood that 



dz d^ 



[z ,x 



(6.31) 



Note the lower limit Xn 



-oo in the above integration, which is equivalent to choosing 



the retarded prescription in 1/9"* 
, , 1 



{X 



\y-) - / 



dp+ e*P^(^" 



-y-) 



2tt p+ + ie 



—i9{x — y 



(6.32) 



This represents our gauge fixing prescription: the LC gauge which we get in this way is 
defined by (cf. eqs. ( |6.25| ) and ( |6.31D ): 



0, 



lim 



0. 



(6.33) 



(This is unique since any ^'-independent gauge transformation would violate the above 
asymptotic condition on a~ .) Equivalently, the ie prescription for the axial pole in the 
LC gauge propagator ( 6.30|) is the retarded prescription shown in eq. (|3.84 ). 

This particular gauge condition is chosen for consistency with the classical solution in 
Sect. 2.3 : There, we have assumed the tree-level source p to have support at x~ > 0, 
and we have chosen the boundary conditions in such a way that the field A^ is located at 
positive x~ as well (cf. eq. (|2.31| )). These assumptions are consistent with the quantum 
evolution provided we use the retarded prescription in the LC gauge propagator of the 
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semi-fast gluons. Indeed, with this choice, the induced charge a comes out with support 
at positive, and relatively large, x~ (see, e.g., eqs. ( ^.16|) -( p.l71 )), which then implies a 
similar property^ for the induced field 5A^ (cf. the discussion after eq. ( |3.72| )). 

To summarize, the LC gauge propagator is obtained as: 



or, more explicitly, 

Qij ^ 



(6.34) 

(6.35) 
(6.36) 
(6.37) 
(6.38) 



where we use the convention that the derivatives written on the right act on the functions 
on their left; e.g., d~Fd~ = d~d~ F{x,y). Also, the operators written on the left 
(right) of G^'^ are regularized with a retarded (advanced) prescription. For instance. 



G-' = 




hi ^ /<++ ^ 'T>)\i 

d+ d+ 


G'- = 






G— = 


-'^'^^ 


d+ ' 



d+ d+ 



dt dt' 



(6.39) 



where l/dj^ has been defined in eq. (|6.32|) , while 



{X 



\y-) 



dp 



+ pip^i^ -y ) 



X 



27r p+ — ie 

The following symmetry property has been used in writing the equations above: 

1 , , , , 1 



\y-) 



id 



\x ) . 



(6.40) 



(6.41) 



R 



Together with the previous expressions for the scalar propagator G, eq. ( |6.13D , and for 
the temporal-gauge propagator G^^ , eq. ( |6.22|) , this concludes our construction of the 
LC-gauge gluon propagator in the background field. 

^"'^If, say, an advanced ie prescription was used instead, the corresponding <j would be located at negative 
x~ , so that a field 5A^ would be generated at a;~ < by quantum evolution, in contradiction with the 
original boundary conditions at tree-level. 
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7 Conclusions 



In this paper, we have given a careful derivation of the JKLW renormahzation group 
equation describing quantum evolution in the high-density regime at small x, and pre- 
pared the framework for the explicit calculation of the coefficients in this equation, to be 
presented in Paper II. 

In this context, we have introduced the physical interpretation of the saturation regime 
at small x as a Color Glass Condensate. This interpretation is supported by the effective 
classical theory generated by the JKLW equation ("the McLerran-Venugopalan model"), 
where (a) the average over quantum fluctuations is replaced by a classical average over 
disorder (as typical for a glass) and (b) saturation is seen explicitly, and involves color 
fields as strong as A* ~ 1/(7 (a value typical for a condensate). 



The derivation presented here differs from the original one in Refs. [pi)| , [Till in that we 
have considered directly the field correlators, rather than the correlators of the fiuctuations 
in the color charge. We thus have explicitly verified that, to the order of interest, the 
quantum correlations at small x are correctly reproduced by the classical MV model 
supplemented with the JKLW evolution equation. 

In our analysis, we have addressed and clarified several subtle points which are essential 
to make the effective theory and the JKLW equation unambiguous. 

To properly formulate the functional RGE, and also the solution to the classical source 
problem, we found it necessary to have an extended structure in x~ . This allows us 
to properly define the functional differentiations with respect to p, as well as various 
expressions arising from quantum fluctuations which appear in the RGE. 

The spatial distribution of the source strength p{x^,Xi_) in x~ has been found to be 
gauge dependent. With the retarded gauge prescription which we have chosen here, we 
avoid initial state interactions of the gluon fleld, and the source has support at positive 
x~ only. In the limit where the source strength is weak, the general RGE reduces to 
the BFKL equation independently of the gauge-flxing prescription, and the longitudinal 
structure is not important. 

To make the classical solution tractable, we have expressed our results in terms of 
a source appropriate for the covariant gauge background fleld problem. When we do 
quantum corrections in this gauge, we found that it was necessary to include a classical 
polarization term in the renormahzation group equation. 

To cope with the temporal non-locality of the action, the quantum version of the 
McLerran-Venugopalan model has been formulated along a Schwinger-Keldysh contour in 
complex time. The structure of this contour has been shown however not to be important 
in the approximations of interest. 

Since the background is independent of time but inhomogeneous in x^, we found it 
convenient to integrate the quantum fluctuations in layers of rather than of . In 
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this way, the longitudinal momenta remain unrestricted, which has allowed us to exploit 
the geometry of the problem and construct the background field propagator. With our 
retarded axial prescription, we have obtained compact and relatively simple expressions 
for this propagator. 

In the sequel to this paper [Q, the expressions derived here will be used to explicitly 
evaluate the coefficients a and x in the RGE, to all orders in the strong background 
field which represents the Color Glass Condensate. We shall also show that, in the large 
Nf. limit, our results become equivalent to the non-linear evolution equations previously 
proposed by Balitsky and Kovchegov |]2l| . 
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